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Abstract 

We explore a model of the world based on real-vector-space quantum theory. In our model the 
familiar complex phase appearing in quantum states is replaced by a single binary object that we 
call the ubit, which is not localized and which can interact with any object in the world. Ordinary 
complex-vector-space quantum theory can be recovered from this model if we simply impose a 
certain restriction on the sets of allowed measurements and transformations (Stucckclbcrg's rule), 
but in this paper we try to obtain the standard theory, or a close approximation to it, without 
invoking such a restriction. We look particularly at the effective theory that applies to a subsystem 
when the ubit is interacting with a much larger environment. In a certain limit it turns out that the 
ubit-environmcnt interaction has the effect of enforcing Stueckelberg's rule automatically, and we 
obtain a one-parameter family of effective theories — modifications of standard quantum theory — 
that all satisfy this rule. The one parameter is the ratio s/uj, where s quantifies the strength of the 
ubit's interaction with the rest of the world and uj is the ubit's rotation rate. We find that when 
this parameter is small but not zero, the effective theory is similar to standard quantum theory 
but is characterized by spontaneous decoherence of isolated systems. 



I. INTRODUCTION: REAL AND COMPLEX QUANTUM THEORY 



Standard quantum theory is based on a complex Hilbert space: density matrices, observ- 
ables and reversible transformations are all represented by linear operators on such a space. 
However, it has been known since the early days of quantum mechanics that many features 
of the theory are shared by two alternative, hypothetical theories, in which the complex 
Hilbert space is replaced by either a real or a quaternionic Hilbert space. For example, in 
their analysis of the logical structure of quantum theory in 1936, Birkhoff and von Neumann 
noted explicitly that their postulates, which were intended to capture this logical structure, 
would be satisfied just as well by the real and quaternionic models as by the standard com- 
plex theory [1]. To be sure, from an empirical point of view there has hardly been any 
contest among these three theories: the complex version has survived every test, and no 
experiment has been done that requires either the real theory or the quaternionic theory 
for its explanation. Nevertheless, over the years researchers have sought a more fundamen- 
tal understanding — not merely an empirical understanding — of the origin of the complex 
structure 

One avenue of investigation along these lines was carried out around 1960 by Stueckelberg 
[3l H]. He began with the representation of probabilities as squares of real amplitudes, 
treating this step as a natural generalization of ordinary probability theory. He felt that 
one needed to explain the complex structure somehow, starting from a real Hilbert space. 
In order to provide such an explanation he imposed the requirement that the theory admit 
an uncertainty principle of a specific form. This requirement led him to introduce a special 
operator J to be used in the expression of the uncertainty principle. The operator J has the 
property that its square is the negative of the identity operator [32j, and the uncertainty 
principle then holds if one requires every observable to commute with J. With this restriction 
on the observables, the theory becomes equivalent to standard complex quantum theory, with 
Stueckelberg's special operator playing the role of the complex number i. We will say that 
a real operator satisfies "Stueckelberg's rule" if it commutes with J. 

In the present paper we take the real-vector-space theory seriously as a potential theory 
of nature and, like Stueckelberg, we consider the possibility that the complex structure is 
somehow to be located or imbedded in the real theory. However, we do not want simply 
to impose Stueckelberg's rule. Rather, we ask whether the complex structure might emerge 
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dynamically in a particular model. To get started we re-express Stueckelberg's construction 
in the setting of basic quantum theory with a finite-dimensional state space. The finite- 
dimensional case requires a step that was not needed in Stueckelberg's work: we have to 
double the Hilbert-space dimension. That is, to model an ordinary quantum system with 
an A'"-dimensional state space, we need a real Hilbcrt space of 2N dimensions. 

Suppose, for example, that one wants to describe only the spin of a spin-1/2 particle. 
A real state vector in a two-dimensional state space is clearly inadequate. To get all the 
allowed quantum states, one needs to double the dimension to four. But a four-dimensional 
real vector space is not the same as a two-dimensional complex space: one needs three real 
numbers to specify a rank-one projection operator (a pure state) in the real space, but only 
two real numbers to specify such an operator in the complex space (these could be the two 
angular coordinates of the Bloch sphere). In order to restrict the set of states, one can 
impose a version of Stueckelberg's rule, namely, that all density matrices commute with the 
matrix 

/ -1 \ 
-1 
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which has the property that {J <S> I)^ — —I <S> I- (We find it convenient to use the symbol 
J for the 2x2 matrix rather than for the larger matrix as in Stueckelberg's papers or in 
the mathematical literature on complex structures.) As we will see in Section II, the theory 
that results from applying this restriction not only to states but also to measurements and 
transformations is equivalent to standard quantum theory for this system, and a similar 
result holds for any other system, regardless of the dimension of the state space. One 
always needs to double the dimension in order for the real-vector-space theory to be able 
to accommodate all the states and operations of the complex theory, and then one needs to 
impose a restriction so as to limit the sets of states and operations in the right way. 

Now, doubling the dimension corresponds to adding to the system a single binary quan- 
tum object. So we can say that the spin of a spin-1/2 particle, that is, a qubit, is equivalent 
to a binary real-vector-space object — a rebit — together with an auxiliary rebit, such that 
the whole system obeys Stueckelberg's rule. One might initially think, then, that to describe 
the spins of n spin-1/2 particles, one would need, in addition to the n basic rebits, another n 
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auxiliary rebits to turn them all into qubits. But it is clear that this is not the case; a single 
auxiliary rebit is sufficient for the whole system, because it is all that is needed to double 
the dimension. This fact has been noted by a number of authors [TBHT^ . For example, it 
has been shown that a quantum computation using n qubits can be simulated efficiently by 
a real- vector-space quantum computation using only n + 1 rebits [T7] . 

We are thus led to consider the following model. Every system is to be described as a 
quantum object in a real vector space, with the same dimension it would normally have in 
the complex theory, and in addition, there is a single auxiliary rebit. We call this auxiliary 
rebit the universal rebit, or ubit, because in this model it needs to be able to interact 
with every object in the world. By invoking Stueckelberg's rule, we could, not surprisingly, 
make our model equivalent to standard quantum theory, as we will show in Section II. 
However, as we have said, in this paper we want to take a different tack: we ask whether we 
can arrive at ordinary quantum theory, or an approximation to ordinary quantum theory, 
without invoking Stueckelberg's rule. If the ubit is interacting with everything, then no local 
observer will be able to control its interactions with distant objects. It is conceivable that 
this uncontrollability could lead to an effective theory that approximates ordinary quantum 
theory, even though the underlying theory is the real-vector-space theory. Whether such an 
approximation is possible is the question with which we begin our investigation. 

One finds that with no further assumptions, a random interaction between the ubit and 
a large environment does not reproduce standard quantum theory. Rather, the ubit quickly 
factors out of the system and becomes irrelevant, and one is left with ordinary real- vector- 
space quantum theory with no ubit, a theory that is in serious conffict with experiment. 
However, the results we present in this paper indicate that if the ubit is rotating sufficiently 
rapidly in its two-dimensional real vector space (rotation is the only internal dynamics 
possible for this simple system) , then one does recover an effective theory that is very much 
like ordinary quantum theory. Our main goal in this paper is to begin to discern the features 
of this effective theory. One feature we might expect to see, and we do indeed see, is that an 
isolated system can undergo spontaneous decoherence with an associated increase in entropy, 
even when no decoherence would be predicted by standard quantum theory. One expects 
such decoherence because, although the system may be isolated in the sense that it does not 
experience any ordinary interactions, the ubit is never isolated and can therefore serve as a 
conduit of information to the rest of the world. 
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It may seem quite fanciful to imagine a special rebit with no particular location, associ- 
ated with the universe as a whole. Indeed in this paper we are not prepared to offer any 
interpretation of this object beyond what the mathematics implies. But we note that ordi- 
nary quantum theory does have a feature that is something like the ubit. For a quantum 
system with definite energy E, even if it is a spatially extended system with many parts, the 
time dependence of the Schrodinger wavef unction is expressed in an overall factor e~*^*/'* 
multiplying the rest of the wavef unction. It is interesting that there is only one such time- 
dependent phase factor for the whole system, not one for each part. (Each part may not 
have a definite energy of its own.) And the phase factor does represent a kind of rotation 
in a two-dimensional real vector space, namely, the complex plane. In our rcal-vcctor-space 
model there is no phase factor, but in its place there is the ubit, with its own internal 
dynamics (that is, the rotation) and its own interactions with other systems. 

One might worry that our model is immediately suspect in that it seems to allow instan- 
taneous communication over an arbitrary distance: a sender Alice could allow her particle 
A to interact with the ubit, which immediately interacts with particle S in a distant galaxy, 
delivering Alice's message to Bob (who has somehow managed to be there). In this paper, 
though, we focus particularly on a limiting case in which it seems that such instantaneous 
communication cannot occur. There are three relevant frequency scales in the model: (i) 
the rotation rate of the ubit, (ii) the typical strength of interaction between the ubit and 
the large environment, and (iii) the typical frequency scale of the local dynamics. We focus 
our attention on the case in which the first two of these frequency scales approach infinity 
with a fixed ratio, while the third remains finite. In this limit we use a heuristic argument to 
obtain a reduced dynamics of the local system. Under our reduced dynamics we automati- 
cally recover Stueckelberg's rule for the states and transformations, which we show prevents 
instantaneous signaling through the ubit. The ratio between the strength s of the ubit's 
interaction with the environment and the ubit's rotation rate ou serves as a single parameter 
that characterizes the effective theory. When that ratio is zero, our results indicate that one 
recovers ordinary quantum theory. Our primary interest in this paper is the case when the 
ratio is small but not zero. 

Though our model does not allow instantaneous communication under the conditions we 
consider, the model itself is manifestly nonlocal and quite contrary to the spirit of special 
relativity (not to mention general relativity). At any given value of the universal time 
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coordinate, the ubit's state changes everywhere in the same way. One can imagine making 
the model local by replacing the ubit with a "ubit field." It is interesting to ask whether 
such a change would ruin the partial agreement we will find with standard quantum theory. 
As wc will sec, it is important in our model that the ubit be interacting with a large 
environment. Over any short time interval, a ubit field at a given location would interact 
with an environment of limited extent, so the "beneficial" environmental effects that we 
rely on would also be limited. We do not explore this question in the present paper but 
confine our attention to the simple model with a single, binary ubit. Despite the underlying 
nonlocality, the fact that we can get a local effective theory makes it seem worth exploring 
the model to see where it leads. 

The paper is organized as follows. In Section 11 we specify what we mean by "real- 
vector-space quantum theory," and we demonstrate the equivalence with ordinary quantum 
theory when Stueckelberg's rule is imposed with no further restrictions. (We do not impose 
Stueckelberg's rule in the later sections.) In Section III we investigate, numerically and 
analytically, the dynamics of the ubit interacting randomly with an environment. We find, 
among other things, that any component of the ubit state that fails to commute with J 
quickly decays. Our next step is to restrict our attention to the limiting case described above 
in which this decay happens instantaneously and continually. In Section IV we consider one 
specific physical example, the spin of a spin- 1/2 particle precessing in a magnetic field, and 
we explore its behavior numerically when the problem is recast in the ubit model. We 
identify three respects in which this behavior deviates from standard quantum theory: (i) 
the frequency of precession is reduced; (ii) there is a long-term dephasing (mentioned above); 
(iii) there is a periodic variation in the purity of the spin state, indicating that the spin is 
periodically becoming correlated and then uncorrelated with the environment. In Section 
V we explain all these effects analytically using perturbation theory, taking the ratio s/uj 
as our small parameter. We find, though, that at least to second order in this parameter 
we can eliminate the strange oscillation in purity simply by reinterpreting the theory. The 
reinterpretation — which does not eliminate either the reduction in precession frequency or 
the decoherence — is presented in Section VI. We consider in that section systems with higher 
state-space dimension than a spin- 1/2 particle, but in this paper we do not analyze the 
higher-dimensional case in detail. Section VII focuses on the fact that in the ubit model 
there is not a unique mapping from the complex theory to the real theory. The choice of 
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mapping amounts to an additional specification of the dynamics beyond what is determined 
by the Hamiltonian. It turns out that one particular choice would render the retardation 
in the evolution unobservable — it would slow all processes by the same factor — leaving only 
the decoherence as a potentially observable effect of the ubit model. We discuss our results 
and draw conclusions in Section VIII . 



II. REAL- VECTOR-SPACE QUANTUM THEORY WITH FINITE DIMENSION 

One can identify four main components of the basic framework of standard quantum 
theory in a complex vector space: (i) States are represented by positive semi-definite op- 
erators with unit trace (density matrices), (ii) A reversible evolution is represented by a 
unitary transformation, (iii) An ideal repeatable measurement is represented by a set of 
orthogonal projection operators Ilj, such the probability of the outcome i when the state 
is p is Tr(njp), and when outcome i occurs, the final state of the system is proportional to 
UipUi. (iv) The state space of a composite system is the tensor product of the state spaces 
of the components. Other kinds of evolution and measurement are certainly possible, but 
they can be derived from the above kind by applying these rules to a larger system and then 
considering the effects on a subsystem. 

The analogous statements for real-vector-space quantum theory are exactly the same, 
except that all the operators are real. In particular this means that a reversible evolu- 
tion is represented by an orthogonal transformation, which is the real version of a unitary 
transformation. 

Let us now write down the differential equation governing the evolution of a state in 
real-vector-space quantum theory. In the usual complex theory, we can write the equation 
of evolution as 

f^^[-zH/h,p]. (2) 

In the real-vector-space theory there is no direct analogue of the Hamiltonian H, but we can 
replace the antihermitian operator —iH/h with an antisymmetric real operator S, so that 
the evolution equation becomes 

fHS.Pl. (3) 
We will call the operator S the "Stueckelbergian" of the system. If the Stueckelbergian is 
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constant, as we will always assume in this paper, then the general solution of Eq. ^ is 

p{t) = e^V(0)e-^*. (4) 

As S is antisymmetric, the operator e"^* is orthogonal. 

We now show how one can recover standard quantum theory from the real-vector-space 
version by adding the ubit and imposing Stueckelberg's rule. Much of what follows in 
this section (minus the interpretation in terms of a ubit) is similar to the account given in 
Ref. m- 

Suppose the system we want to describe has a d- dimensional (complex) Hilbert space. 
Then we start by considering a (i-dimensional real-vector-space object A together with the 
ubit U. Consider any matrix M that might apply to the UA system, whether it be a 
density matrix, an orthogonal evolution operator or antisymmetric Stueckelbergian, or the 
projection operator associated with a measurement outcome. The 2d x 2d matrix M can be 
written as 

mJ"'""'A, (5) 
\^ Mio Mn J 

where each Mj-^ is a real d x d matrix and the subscripts refer to a pair of orthogonal basis 
vectors in the ubit's two-dimensional space. We now impose Stueckelberg's rule: we insist 
that M commute with Ju®Iai where Ju is the 2x2 matrix we mentioned in the introduction 
and I A is the dx d identity. (We include the alphabetic subscripts to indicate which system 
the operator acts on.) That is, we insist that 

Moo Moi W -Ia\ I -Ia\I Moo Moi 
Mio Mn y \ /a y 1^ /a y \^ Mio Mn 
which implies that Moo = Mn and Mio = —Moi. Thus we can write 



(6) 



, Moo -Mio . 
M = °° '° . (7) 

Mio Moo 

We can map any such matrix into a smaller, dx d complex matrix, such that under this 
mapping, the laws of real-vector-space quantum theory become the laws of complex-vector- 
space quantum theory. The mapping is this: For a matrix M representing an orthogonal 
transformation, a Stueckelbergian, or a projection operator, we have 

, Moo -Mio , 

M = I I ^ >1 = Moo + ^Mio. (8) 

Mio Moo 
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And for a density matrix p, we have 

a = 2{poo + ipio). (9) 



Poo — Pio 
Pio Poo 



The special rule for density matrices is simply to make sure every real or complex density 
matrix has unit trace. Note that pio must be an antisymmetric matrix in order that p be 
symmetric. This implies that Tr p = Tra, since pio is traceless. 

One can verify that Eq. ([s]) faithfully preserves matrix multiplication: If Mi Aii and 
Ma M2, then M1M2 M1M2. And Eqs. g and Q together preserve the trace of 
a density matrix times a projection operator: If p — > cr and 11 — T, then Trllp = TrTa. 
These two facts guarantee that the real- vector-space laws (restricted by Stueckelberg's rule) 
are equivalent to the complex-vector-space laws under this mapping. 

So far we have been imagining going from the real theory to the complex theory. One 
can just as easily go the other way around. For example, given a complex Hamiltonian if, 
we can write the corresponding Stueckelbergian 5* as 

( Rei-iH/h) -lm(-iH/h)\ 
S = \ ^ ' ' ^ ' M =Iu®Re{-iH/h) + Ju®\riy{-iH/h). (10) 
\ lm{-iH/h) Re{-iH/h) 1 

Given a complex density matrix a, one obtains the corresponding real density matrix p by 
performing a similar operation [33] : 

I I Re a -Ima \ 1 
p=- \ = -{Iu(»Recr + Ju^lm(r). (11) 

^ \ Ima Re a I ^ 
Note that in this real-vector-space setting, under Stueckelberg's rule, there is no such thing 



as a state vector. In fact the purity Tr p^ cannot be greater than 1/2: from Eq. (11), we 
have Trp2 = (l/2)Tr[(Re a)^ - (Ima)^] = (l/2)Tr(T2 < 1/2. Thus every state has to be 
represented by a density matrix, even if it corresponds to a pure state in standard quantum 
theory, and the minimum rank of any density matrix is two. 

Of course Stueckelberg's rule is crucial here. Without this restriction, one could indeed 
have a state vector in the real-vector-space theory — in a sense such a state would be purer 
than any pure state in standard quantum theory. Also note that a general orthogonal matrix 
in 2d dimensions can be characterized by 2(P' — d real parameters, whereas a unitary matrix in 
d dimensions requires only d'^ real parameters. So hardly any of those orthogonal matrices 



correspond to unitary matrices. This is one sense in which the unrestricted real-vector- 
space theory allows too many possibilities. In the following sections we will not impose 
Stueckelberg's rule but will try to achieve its effects in another way. 

It is useful to note that, in the absence of any restrictions, every real matrix M acting 
on the UA system can be broken uniquely into two parts Mc and Ma, which respectively 
commute and anticommute with Ju ® Ia- We can write the two parts as 

Me = ^ [M - {Ju ® Ia)M{Ju ® I a)] (12) 

and 

M„ = ^ [M + {Ju ® Ia)M{Ju ® I a)] . (13) 

In these terms one can see, for example, that if all observables and transformations commute 
with Ju ® Ia-i then the anticommuting part of a density matrix, p^, can have no observable 
effects. Let an initial density matrix p be transformed by an orthogonal transformation O 
and then tested for a property represented by the projection operator 11. If O and 11 commute 
with Ju <S> Ia, then the contribution from pa to the probability of the "yes" outcome is 

Tr{UOpaO^) = (l/2)Tr {UO [p + {Ju ® Ia)p{Ju ® Ia)] O^} = 0. (14) 

(In the second term inside the trace, we can move one factor of Ju®Ia through O, 11, and 
so that it combines with the other factor of Ju ® I a to yield —lu ® I a- The two terms then 
cancel.) Thus if we impose Stueckelberg's rule on all observables and transformations, the 
physical predictions of the theory will not depend on whether we also impose this restriction 
on the set of allowed states. 

III. THE UBIT INTERACTING WITH THE ENVIRONMENT 
A. Specification of the model 

Ultimately we want to consider a system A interacting with the ubit U , which is also 
interacting with an environment E (but there will be no direct interaction between A and 
the environment) . We will take the Stueckelbergian of the entire system to be of the form 

S = -ujIe ®Ju®Ia + sBeu ®Ia + Ie® Sua, (15) 
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where the subscripts again indicate the system on which each operator acts and we use a 
hat (rather than the subscript EUA) to distinguish those operators that act on the entire 
system. The operator Ju generates rotations of the ubit, so u is the ubit's rotation rate. The 
operator Beu characterizes both the interaction of the ubit with the environment and the 
internal dynamics of the environment itself, and s determines the scale of these interactions. 
Finally the operator Sua is the local Stueckelbergian, the one part of S that we imagine can 
be controlled by an observer. 

This last part deserves some discussion. The ubit is not localized, but we are assuming 
that it is available to be manipulated and measured by any observer. That is, our local 
observer — whom we will call Alice — can arrange for the implementation of an arbitrary 
Stueckelbergian Sua involving the ubit and the local real-vector-space object A and can 
measure any observable on the UA system. One might worry that different physical systems 
all over the universe will be competing to achieve contradictory effects on the same ubit. 
Indeed we will see that something along these lines does happen. It will turn out that the 
interaction of the ubit with the environment severely limits what Alice will actually be able 
to do. However, we do not impose any such restriction in the basic model. 

We now describe our method of generating the matrix Beu- One could reasonably model 
the environment as a collection of, say, rebits or higher-dimensional objects, each having 
some random interaction with the ubit but no interaction with each other. However, when 
we do our numerical experiments, we would like each run to be reasonably reproducible, as it 
would be for a very large environment; so we want Beu to include as many randomly chosen 
parameters as possible without having to make the environment's dimension intractably 
large. We therefore model the environment as a single system, with the matrix Beu simply 
chosen at random. More precisely, taking to be the Hilbert space dimension of the 
environment, we create a 2N x 2N matrix R, each component chosen uniformly between 



— 1 and +1, which we then antisymmetrize and multiply by ^JQ/N . The factor l/\fN 
guarantees that the typical size of an eigenvalue of Beu will not depend on A^. The factor 
v^G has been inserted for later convenience. Thus Beu = i^/Q/N)iR - R^)/2. We assume 
that this matrix is written in a tensor-product basis of the environment and the ubit, with 
the ubit basis being such that the Ju of Eq. (15) has the standard form given in Eq. ([T]). 
(A mere rotation of the ubit basis does not change Ju-, but Ju would pick up a minus sign 
under a reflection.) For most of our numerical runs the environment dimension A^ is 200. 



11 



To get some insight into how the EUA system will evolve, in the following two subsections 
we restrict our attention to the simple case in which the system A has no dynamics of its 
own and is not interacting with the ubit. That is, the local Stueckelbergian Sua is zero, 

so that S = Seu ® Iai where Seu = —^Ie ® Ju + sBeu- We also assume that Alice has 
prepared the UA system in an initial state that is uncorrelated with the environment, so 
that the initial state of the whole system is of the form 

p(0) = pE® PUA- (16) 

The state puA, which could be pure because we are not imposing Stueckelberg's rule, and 
which may exhibit entanglement between the ubit and the A system, can always be written 
as 

PuA ^lu® a^^^ + Ju® a^^^ + Xu® a^^^ + Zu ® a^^\ (17) 

where Xu and Zjj are the usual Pauli matrices acting on the ubit's space and the a's are 
operators on As space. With Sua being zero, the matrices {a^^\ a'^'^^ a^^\ a'^^''} will not 
change. That is, at a later time t, when the state of the whole system is p(t), the state of 
the U A system will be of the form 

PuAit) = T^Epit) = Iu® a(^) + u^'\t) ® a(^) + mW(^) ® "^""^ + ^^^H^) ® a^^^- (18) 

(The part will not change, since lu commutes with the Stueckelbergian.) Our aim here 
is to follow the evolution of the ubit matrices u^'^\u^'^\ and u^^\ 

The evolution of these matrices will depend to some extent on the initial state Pe of 
the environment. Numerically we have tried both a randomly chosen pure state and the 
completely mixed state, and we have found that for a sufficiently large dimension of the 
environment's Hilbert space, the results are almost the same in both cases though they tend 
to be somewhat smoother in the latter case. For simplicity, then, in all of our calculations we 
will choose the the initial environment state to be the completely mixed state pe = Ie/N . 

Under these assumptions, our numerical results can be summarized as follows. The 
function u'^'^^ is of the form ^{t)Ju, where 7(t) is an initially oscillating function that finally 
approaches a constant value between zero and one. Thus the Ju part of the state diminishes 
but does not disappear. On the other hand, the functions u^-^^ and u^^^ both become linear 
combinations of Xu and Zu whose coefficient vectors rotate in the X-Z plane and finally 
decay to zero (at least to a very good approximation, which we expect will be exact in 
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the limit of an infinite-dimensional environment). The rotation of the coefficient vectors is 
simply a manifestation of the rotation of the ubit in its two-dimensional real state space. 
That M^^) and u'^^^ decay to zero tells us that the state of the U A system eventually obeys 
Stueckelberg's rule, since the remaining operators lu and Ju commute with Jj/. These 
results can be understood through perturbation theory, as we now show. 



B. The function j{t) 

The 2x2 matrix u^"^^ can be written as 

u^'^\t) = ^Tr^ [e^^^*(/£ ® J[/)e-^^^*] . (19) 

As we have said, this matrix remains proportional to Jjj — this follows from the fact that J 
is antisymmetric while /, X and Z are symmetric — so that u^'^\t) — ^{t)Ju for some real 
function ^{t). Here we try to estimate ^{t). We can write it as 

l{t) = -^Tr [e^^^*(/E ® Ju)e-'^''\lE ® Ju)] ■ (20) 

To apply perturbation theory, it is helpful to define a Hermitian matrix G as follows: 

G = — - -i/fi ® Jt/ + -{iBeu). (21) 

U! CO 

We can think of G as 

G = Go + XV, (22) 

where 

Go = —He <8) Ju and V — IBeu, (23) 
and X — s/uiis our perturbation parameter. In terms of G, we have 

7(^) = -^Tr [e-'^'^'ilE ® Ju)e'^''\lE ® Ju)] ■ (24) 

Note that Go, which is a 2N x 2N matrix, has only two eigenvalues, +1 and —1, each 
associated with an A'"-dimensional subspace. Let P+ and P_ be the projection operators on 
the subspaces corresponding to the eigenvalues and —1 respectively. We can write these 
operators explicitly as 

P+ = ]^Ie ® {h - iJu) and P_ = ^/^ (g) {h + iJu) ■ (25) 
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To do degenerate perturbation theory, we choose a basis that diagonahzes the matrix V in 
each of these two subspaces. Let |$+) and |$~) be the elements of this basis. That is, the 
vectors |$^) are the eigenvectors of = P_|_VP+, and the vectors 1$") are the eigenvectors 
of = P^V Thus ($^1^1$^) is zero if m 7^ n but ($^1^1$^) is typically nonzero. 

Before proceeding, it is worth noting certain symmetries that follow from the fact that 
any Stueckelbergian is a real, antisymmetric matrix and that P_ is the complex conjugate of 
P+. First, the Hermitian matrix V~ is the negative complex conjugate of . From this it 
follows that we can take I'l',^) to be the complex conjugate of |$^), and if f„ is the eigenvalue 
of associated with |$^) then —Vn is the eigenvalue of V~ associated with |$~). Similarly, 
the eigenvectors of the Hermitian matrix G can be written as I^E'^) and |^E'^), corresponding 
to eigenvalues Qn and —Qn respectively, where |\1'^) is the complex conjugate of \^n)- 



In the expression (24) for 7(t), we make the substitution 

N 



n=l 



Leaving the eigenvalues unanalyzed for now, we use standard time-independent perturbation 
theory to write the eigenvectors in terms of the unperturbed eigenvectors 1$^) and |$~). 
Specifically, we expand and |\E'~) to second order in A (because there is no first-order 



contribution to 7(t)) and insert this expansion into Eq. (26), which in turn is inserted 
into Eq. (24). The second-order expansion of |^E'^) is given in Appendix A. The resulting 



expression for 7(t) comes out to be 



.2 N 



7(t) = l-^ E m\V\K.)?{^'^o,[{gn + g^)ujt\). (27) 

n,m=l 

With a couple of assumptions about the matrix elements ($^|\^|$~) and the distribu- 
tion of values of gn-, we can obtain an explicit functional form for 7(t). First we assume 
that because many random values are being summed to get the squared matrix element 
we can replace this factor, for each value of m and with its ensemble aver- 
age, that is, the average over all possible matrices Beu generated by the random procedure 
specified above. To find this ensemble average, we begin with 

1 ^ 1 

K\y\^-J?) = ^, E (K*nim;;)r> = ^(Tr(p+i^p-i^)), (28) 

n,m=l 



where the angular bracket indicates the ensemble average. Inserting the expressions (25) into 



Eq. (28) and writing V in terms of the random matrix P, we find that each term involving 
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Ju is zero because it is proportional to the ensemble average of the product of two distinct 
elements of the random matrix R. The remaining terms give us 



($+|l^|$-)r> 



1 



4Ar2 



(Tr(F2)> 



16iV3 



(29) 



where we have used the fact that the average square of a component of R is 1/3, and we 
have neglected (Tr(i?^)) and (Tr[(/2"^)^]) because they are of order whereas (Ti{R^R)) is 



of order A^^. We will use the value 1/A^ in place of |($+|V"|$„)p in Eq. (27). 



We now turn to the eigenvalues ±(7^ of G. Recall that the unperturbed eigenvalues, that 
is, the eigenvalues of Go? are simply +1 and —1. To lowest nontrivial order in A, we can 
write Qn = ^ + Xvn- (Again, t>„ = is an eigenvalue of V'^.) Because comes 

from linearly transforming a random matrix, for large A^ we expect its eigenvalues to 
approximately exhibit a semicircular distribution. To write this distribution explicitly, we 
find the ensemble average of (1/A^) X] "^n^ reasoning as above but now with a smaller matrix: 



N 



N 



Tr 



AN 



Tr {V^)) ^ 1. 



(30) 



Let ri{v)dv be the expected number of eigenvalues of V'^ lying between v and v + dv, so 
that the normalization of r) is fixed by the condition p'"""" ri{v)dv = N, where Vmax is the 
maximum value of v. Then the unique semicircular distribution satisfying (f^) = 1 is given 
by 

^{v) = {NIt^)^\-{vI2)\ (31) 



so that Vr, 



2. To get our analytic expression for 7(t), we replace gn in Eq. (27) with 



1 + Af and gm with 1 + Af', and we assume both v and v' are distributed according to 77. 
With these approximations, we have 



7(t) = 1 - A^ + 



A^ 

ATS 



ri{v)ri{v') cos [(2 + Xv + Xv') ut] dv dv' . 



(32) 



The integral can be done exactly, and we obtain 



lit) 



X' + 



Ji{2st) 



(33) 



where Ji is the Bessel function of order 1. When t is large, 7(t) approaches the constant 
value 1 — A^. That is, the J part of the matrix pjjA is reduced but does not vanish. This 
analytic expression is compared with our numerical results in Fig. [1} 
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FIG. 1: The evolution of the coefficient 7(t) of Ju under the influence of the environment. The 
initial value of 7 is set to unity, and the system is initially in a product state, with no correlation 
between the ubit and the environment. The dots (blue) show the numerical results, and the curve 



(red) shows our result from perturbation theory (Eq. (33)). Here = 200, s = 10 and a; = 100. 
(Time is measured in arbitrary units. If t is interpreted as being in zs, for example, then s = 10 
zs~^ and oj = 100 zs~^.) In this figure, as in every figure in this paper, the numerical results are 
for a single run of the simulation. The results from one run to the next are extremely consistent. 

C. The matrix u^^) 

We would now like to get an analytic expression for u'^^\ which starts out as X at t = 0. 
The matrix u^^^ can be written as 

u^^\t) = ^Tte [e^^^*(j£; ® Xu)e-^^^'] . (34) 

We can always write a 2 x 2 matrix as a linear combination of basis elements, so we can 
write u^^'^ as 

«W(t) = /3(^)(t)/ + l3^J\t)J + (3^''\t)X + (3^^\t)Z, (35) 

where the /3's are real-valued functions. As X is symmetric, (3'^'^^ must equal zero since J is 
anti-symmetric. We also know that equals zero since the identity commutes with e^'^^^: 

P^'^ = ^Tr [{Ie ® Iu)e'^''\lE ® Xf;)e-^-^*] = ^Tr [Ie ® Xu] = 0. (36) 
Thus u^^\t) = P^^\t)X + /3(^)(t)Z. 
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To calculate (3^^^ and f3^^\ we again apply perturbation theory. We begin with the 
following expression for (3^^^: 

P^'^Kt) = ^Tr [{Ie ® Xt;)e^--*(/s ® Xu)e-'^^'] . (37) 

When writing e"^-^*^* in terms of the unperturbed eigenvectors, it is convenient to factor 
each unperturbed eigenvector into a tensor product of the environment and ubit parts. The 
eigenvectors of —iJ (the ubit part of Gq) are |+) = and |— ) = -^{\), corresponding 

to the eigenvalues +1 and —1. Thus, 

l-^-^) = 10+) ® 1+) (38) 

\K) = \<l>n) ^ H (39) 

for some environment states (where is the complex conjugate of |0^)). We now 
write X in terms of |+) and |— ): 

X = 




i(l+)(- 



(40) 



Inserting Eqs. (26) and (40) into Eq. (37), we find an expression for (3^^^ to zeroth order in 



A for the eigenstates. There are no first order terms, and the zeroth-order terms turn out 
to be sufficient to give us good agreement with the numerical calculations. The resulting 
expression comes out to be 

TV 

/3W(t) = 2 cos [{g^ + g Jut]. (41) 

n,m=l 

Using the same assumption as before about the distribution of the values of gn, we can 
find an expression for j3^^^ independent of the details of the eigenstates and eigenvalues 
of G. We also assume that for a sufficiently large environment, we can approximate each 



I (0+ 1 by its ensemble average Plugging this value back into Eq. (41) and assuming 



the same semicircle distribution as before, we arrive at the final expression: 



Ji(2st) 

St 



1 2 



cos(2c(;t). 



(42) 



In a similar way, we find that 



1 ^ 

f^^'^^^^ = N ^ M\<PrJ\'sm[{gr, + gr^)u;t] 



n,m=l 



Ji{2st) 



n 2 



St 



sin(2u;t). 



(43) 
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The functions (3^^'' and together show what happens as u^^"* evolves over time. The X 
and Z parts of the ubit's matrix u^^^ rotate into each other and eventually decay to zero. 
Fig. M compares our analytic expressions for both and (3^^^ to our numerical results. 

Solving for u^^'^ yields a similar result. The ubit matrix u^^^ begins as Z, and the Z and 
X parts again rotate into each other as they decay down to zero. 



coefficient of X 
1.0 




coefficient of Z 
1.0 



time 0.0 




time 



FIG. 2: Evolution of the matrix u^^\t) = + /3(^)Z, which begins as X. Here we plot (3^^^ 

and (3^^^ as functions of time. One can see that the vector {f3^^\ P^^'^) rotates in the plane and 
finally decays to zero. The dots (blue) show the numerical results, and the curves (red) show the 



result of our perturbation theory calculation (Eqs. (42) and (43)). Here = 200, s = 10 and 
CO = 100. 



D. Projecting onto the space of matrices that commute with Seu 

The above calculations show that 7(t), (3^^\t), and (3^^\t) all approach asymptotic val- 
ues. We could have obtained these asymptotic values more directly in the following way. 
First, we can break p(0) into two parts: the part pn lying in the space of matrices that com- 
mute with Seu ® I a, and the part p± perpendicular to this space [31]. Then the evolution 
(again in the absence of any local Stueckelbergian Sua) becomes 

p(t) = e(^^^®^-^)*(p|| + p^)e-('^EU^'A)t ^ + ^iSEU^lA)tp^^-(SEU^lA)t_ (44) 

The second term is what provides the oscillations we saw in the graphs in the preceding 
subsections. Asymptotically the effect of this second term on the UA system disappears as 
the oscillating components tend to cancel each other out when we trace over the environment. 
It is only py that contributes to the final state of UA; so in order to find that state we could 

18 



have focused only on p\\ (to which we could have applied perturbation theory as in the above 
calculations). 

For the remainder of this paper, in our analytical work we will adopt the following ansatz. 
First, we assume that both s and w, that is, the scales of the two terms in the EU part 
of the Stueckelbergian, are much larger in magnitude than the spread in eigenvalues of the 
local Stueckelbergian Sua- We imagine taking the limit as s and uj both go to infinity while 
their ratio remains fixed. As s and uj get larger, the oscillations and decay we observed in 
the preceding subsections simply proceed at a faster rate without changing their form in 
any other way. Consider, then, the evolution of p over any short interval of time. By the 
time Sua has had any appreciable effect, the asymptotic value of puA due to the action of 
Seu "will have already been reached. Therefore, for the purpose of computing puA we will 
(i) ignore any initial p_L and (ii) assume that p||, as it evolves, is continually projected into 
the space of matrices that commute with Seu ® I a- That is, we will use this continual 
projection in place of the exact evolution due to Seu in all our later analytical calculations. 
(But in our numerical work we will follow the exact dynamics.) Note that this method 
of continual projection does not necessarily provide a good approximation to p itself. The 
rapidly oscillating part p_L does not go away, but it is irrelevant for computing puA- To 
remind ourselves that we are dealing only with py and not with the full density matrix 
p, we will keep the subscript "||" when referring to the state of the whole system. In 
the following paragraphs we write down some general consequences of the assumption of 
continual projection. 

First we modify the equation of evolution, Eq. ([3]), so that it does not allow p\\ to evolve 
away from the space of matrices that commute with Seu ® I a- The modified equation is 

'^-^ = V{[lE®SuA.V{p\0\), (45) 

where V projects onto this space. (The second V is unnecessary, since p\\ is already in the 
space into which V projects. We include it only so that the equation would preserve the 
trace of any density matrix.) Assuming that Seu is non-degenerate, the only matrices that 
commute with it are linear combinations of projections onto its eigenstates. Thus we can 
express the action of "P on a generic matrix M as follows: 

V{M) = J2 ® T^^Eu m,){^,\ ® Ia) M] . (46) 

j 
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Here we use the index j to stand for the combination of n and ± in (Again, the 

vectors are the eigenvectors of Seu-) The most general form of py as a function of time 
(still assuming that Seu is nondegenerate) is 

j 



where (Jj{t) is a matrix acting on the space of the A system. Inserting Eqs. (46) and (47) 



into Eq. (45), we get the following equation for the evolution of the cr/s. 



^ = [{^j^\I^®SuA\^j),a,], (48) 

where the quantity {^jIIe ® SuaI'^j) is to be interpreted as a "partial expectation value," 
in which the EU vector |\l'j) combines with the EU part of Ie ® Sua to leave a matrix that 
acts on the space of the A system. 

Now, the local Stueckelbergian Sua must be a linear combination of the four ubit matrices 



/[/, Ju, Xu, and Zu, each in a tensor product with some matrix of the A system. Eq. (48) 
thus calls on us to evaluate the quantities {'^j\lE^Iu\'^j), i'^jl^E^Jul'^j), {^jUs^Xul'^j), 
and {^j\lE ® Zul'^j). The first of these is clearly equal to unity. If |\E'j) were simply a 
random state, then the other three would have typical values that diminish in magnitude 
proportional to I/a/JV for large N [22] • But |\E'j) is not a random state. Rather, it is an 
eigenstate of Seu = —^Ie ® Ju + sBeu where Beu is random. The presence of uIe ® Ju 
in this matrix prevents {"^jIIe ® Jul'^j) from going to zero as approaches infinity, but it 
does not similarly protect {^j\lE ® Xul'^j) or {"^jIIe ^ Zul'^j). We find numerically that 
even with a nonzero value of u these last two quantities have typical values that approach 
zero as I/a/ZV, while the typical size of {^jIIe ® Ju\^j) approaches a nonzero constant. Let 



us define 



z/. 



-t{^j\lE(S)Ju\'^j), (49) 



in which the i has been inserted to make z/,- real. In Section V we will estimate uj, but for 



now we simply use it to rewrite our basic equation (48). According to what we have just 
said, in the large A^ limit we can ignore any part of Sua that is proportional to Xu or Zu, 
so that in effect the most general Sua has the form 

Sua = Iu®La- Ju® Ka, (50) 
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where La is an antisymmetric real matrix and Ka is a symmetric real matrix. Inserting this 



form into Eq. (48), we get 

^ = -i[vjKA + iLa, (Tj\. (51) 
(The matrix aj can be complex, as long as the imaginary parts cancel out when we do the 



sum in Eq. (47).) Evidently the Hermitian matrix VjKA + iLa is playing a role like that 
of H/h, except that because of the j dependence in z/j, different components aj can have 
different effective Hamiltonians. We will get a sense of what consequences this fact has as 
we consider in Sections IV and V the special case of a precessing spin. 

The fact that j\lE ® Xjj\^ j) and {^j\lE ® Zu\^ j) become zero in the large limit 
has another important consequence. First, it means that when we expand as a 

linear combination of the ubit matrices /[/, Jjj, Xjj, and Zu, the environment matrices 
multiplying Xu and Zjj must have zero trace. Therefore, for any py of the form given in 



Eq. (47), the density matrix puA resulting from tracing over the environment cannot include 
any term proportional to Xu or Zu (in the limit as approaches infinity), so that puA 
commutes with Jjj. Thus both our local Stueckelbergian Sjja and our local density matrix 
PuA commute with J[/, and in this sense we have recovered Stueckelberg's rule through the 
interaction of the ubit with the large environment. We have not explicitly ruled out the 
possibility of a measurement operator that does not commute with Ju, but if Alice were 
to manage to perform such a measurement, represented by a projection operator Pua, the 
anticommuting part {PuA)a = {^/'^)[Pua + {Ju^Ia)Pua{Ju®Ia)] would make no observable 
difference because TT[pi/A{PuA)a\ is equal to zero for any puA that commutes with Ju ® Ia- 
We hasten to add, though, that this effective enforcement of Stueckelberg's rule by our 
projection hypothesis does not make our theory equivalent to standard quantum mechanics. 
It eliminates unwanted states and unwanted terms in the Stueckelbergian, but as we will see, 
it is not equivalent to simply imposing Stueckelberg's rule as in Section II. The interaction 
with the environment yields a different effective dynamics. Our main task in the rest of this 
paper is to characterize the differences. 

Finally, one might wonder whether it really makes sense to assume, as we have done, 
that the parameter s, which is the size of a typical eigenvalue of sBeu, is much larger in 
magnitude than the spread in eigenvalues of the local Stueckelbergian. After all, the idea 
underlying our model is that the ubit's interaction with each component of the rest of the 
world should be similar to its interaction with the local system. 

21 



In fact there is no contradiction here. In a more realistic model of the environment, the 
size of a typical eigenvalue of the ubit-environment Stueckelbergian would reflect not just the 
strength of interaction between the ubit and a single component of the environment. It would 
also reflect the size of the environment. Suppose, for example, that the environment consists 
of n rebits and that the ubit interacts with each one via a simple 4x4 Stueckelbergian matrix 
with eigenvalues {i^, i^, — i^, — i^}. Then even if those individual rebits do not interact with 
each other, the square root of the average squared magnitude of an eigenvalue of the whole 
interaction Stueckelbergian is equal to iy/n. Thus the typical size of an eigenvalue grows 
with the size of the environment. This model of the environment as composed of independent 
systems is not the one we have chosen for our numerical simulations, but this argument shows 
that it is reasonable to assume that the matrix sBeu has large eigenvalues: its eigenvalues 
are large by virtue of the large size of the environment, even if the strength of interaction 
between the ubit and any small component of the environment is of limited magnitude. 



E. No signaling 

To summarize the last subsection: by imagining both s and u going to infinity with a 
fixed ratio, we were led to consider only the part of the global density matrix that commutes 
with Seu ® I A, and we assumed that during the evolution, this part is continually projected 
onto the space of matrices that commute with Seu ® Ia- This assumption led us to the 



form (47) of the density matrix, which evolves according to Eq. (48). Next, we considered 
the implications of the environment's dimension becoming infinitely large. (In our model, 
this means that the number of independently chosen random parameters becomes infinitely 
large.) In this limit, we concluded — admittedly on the basis of numerical evidence — that 
certain terms of the local Stueckelbergian will become inconsequential, because the random 



nature of the matrix Beu causes the contributions of these terms to Eq. (48) to vanish. The 



only terms in Seu that can have any effects, then, are those that commute with Jjj. In that 



case Eq. (48) can be written in the specific form (51). We now show that this form of the 
equation does not allow signaling between two observers Alice and Bob if the systems they 
hold, A and B, have no direct or indirect interaction between them except through the ubit. 

We first have to write down what it means that A and B are not interacting except 
through the ubit. In standard quantum theory, two isolated and therefore non-interacting 



22 



systems A and B have a Hamiltonian of the form 

Hab = Ha®Ib + Ia® Hb. (52) 

(We use script letters to refer to complex- vector-space systems.) Converting this Hamilto- 
nian to real-vector-space language as in Section II, we have that the Stueckelbergian is 

hSuAB = /[/ ® Re {-iHa ®Ib- iIa ® Hb) + J[/ ® Im {-iHa ® h - Ha ® Hb) . (53) 



We can write this operator in a form like that of Eq. (50): 



SuAB = (lu ®La- Ju® Ka) ®Ib + Ia® (Iu ®Lb- Ju® Kb). (54) 



Given that we are ruling out any terms proportional to Xjj or Zjj^ the form given in Eq. (54) 
is the most general form possible for a pair of isolated systems (that is, isolated except for 
their interaction with the ubit). 

For the pair AB^ we can rewrite Eq. (E7| as 



where r,- is an operator on the space of the AB system, and the hat here labels an operator 



on the whole EUAB system. With a Stueckelbergian of the form (54), the equation of 



evolution for Tj is a modified form of Eq. (51): 



dr- 

= -i[{ujKA + iLa) ®Ib + Ia® {i^jKB + Hb), tj]. (56) 

We are now in familiar territory. The matrix Tj acts on the space of A and B, but it evolves 
according to an effective Hamiltonian that includes no interaction between the two systems. 
Therefore the partial trace of tj over either of two systems evolves under its own effective 
Hamiltonian, with no influence from the other system. E.g., 

^Trs Tj = -i [ujKA + iLa, Tr^ Tj] . (57) 

It follows that puA, which is 

PUA = Tr^BPii = ^ E^r^ (l*i)(^.l) TrB (r,(t)) , (58) 

j 

evolves independently of what happens to system B. That is, Bob's choice of local Stueck- 
elbergian Sub = Iu ® Lb — Ju ® Kb cannot affect what Alice sees. Now, Bob can perform 
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operations other than simply applying a Stueckelbergian to the U B system for a period of 
time. He can allow system B to interact with other systems (which by assumption are not 
interacting with ^4) and in particular he can make measurements. But any such operation 
can be accounted for simply by expanding the definition of system B. (Bob may observe a 
definite outcome of a measurement, but he is not allowed to convey to Alice any information 
about this outcome. So to describe the state Alice observes we need to keep all the out- 
comes, and we can do this by letting B become entangled with the measuring device with 
no collapse.) We conclude that in this model, in the limiting case we are considering, there 
can be no signaling through the ubit. 

To be sure, if s, c^, and remain finite, then there can be signaling through the ubit. 
It would be interesting to determine quantitatively how the degree of signaling (suitably 
defined) depends on s, a;, and A", but we leave this question for future work. Here we focus 
on the limiting case. 

Thus for any fixed value of the ratio s/u — A, we should get an effective theory that 
is a no-signaling theory. The effective theory will typically not be the same as quantum 
mechanics — it will be a modification of quantum mechanics. The results of the next two 
sections indicate that for the special example we consider there — a precessing spin — we 
recover standard quantum theory when A goes to zero. In Appendix B we present a related 
argument that does not depend on taking either s or A" to be infinite and that applies to a 
general system. There we show that with fixed s and A^, the evolution operator e^* becomes 
equivalent to the standard quantum mechanical operator when ui goes to infinity. 

IV. A PRECESSING SPIN— NUMERICAL SIMULATIONS 
A. The Stueckelbergian and the initial state 

In the numerical work of the preceding section wc considered the interaction of the ubit 
with the environment, in the absence of any interaction with the local system A. We now 
add this interaction and study the simplest possible case, a precessing qubit. For definiteness 
we take the qubit to be the spin of a spin- 1/2 particle in the presence of a constant magnetic 

— * — * 

field B. With B in the positive z direction and the particle having a negative charge, the 
usual Hamiltonian can be written a.s H = h{Q/2)Z, where Q is the angular precession 
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frequency, equal to B times the magnitude of the particle's gyromagnetic ratio, and Z is 
again the Pauli matrix for the z axis. The precession will be in the right-hand sense around 
the positive z axis. 

We can use the correspondence given in Section II to re-express the same phenomenon 



in terms of a rebit A and the ubit. The Stueckelbergian is obtained from H as in Eq. (10): 



Sua = Ju® lm{-iH/h) = -^Ju®Za. (59) 
Let the initial spin state be in the x direction; that is, the initial density matrix of the qubit 



is p = (l/2)(/ + X). To re-express this state in the ubit model, we use Eq. (11): 



Pua{^) = \iu®{Ia + Xa). (60) 

If we were to let the UA system evolve under the Stueckelbergian Sua from the starting 
state pua{^), it would exhibit the standard precession simply re-expressed in real- vector- 
space terms. That is, we would have 

PuA{t) = ^[Iu<®Ia + cos{nt)Iu ®Xa + sm{nt)Ju ® Ja] ■ (61) 

The combination Jjj ® Ja takes the place of the Pauli matrix Y. 

But we are interested in the evolution of the state of the whole system, with the initial 
state 

Pll (0) = -^Ie ® /f/ ® (I A + Xa) , (62) 
under the full Stueckelbergian 

S = -luIe ® Ju ® Ia + sBeu ® Ia + Ie® Sua- (63) 
(As in Section III, we are taking the initial state of the environment to be the completely 



mixed state. We make this choice as a matter of simplicity. Note that Eq. (62) is a special 



case of Eq. (47).) With this initial state and Stueckelbergian, we find the state of the UA 



system at time t by computing 

p,;^(t) = Tr^(e^*/)||(0)e-^*). (64) 

The results of Section III lead us to expect that any component of puAif) involving Xu or Zjj 
will be extremely small, and indeed this is what we observe numerically — those components 
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seem to approach zero as increases. This leaves four basic symmetric matrices in which 
we can expand puA{t)- The expansion can be written as 

PuAit) = ^[Iu®Ia + bJu ^Xa + hyjjj ®Ja + hJu ® Za] • (65) 

Thus in our model we can still imagine the evolution of the spin as a path through the Bloch 
sphere, with the Bloch vector defined as 6 = {h^, by, hz). This vector must have a length no 
greater than unity since puAit) is positive semi-definite. 

As we mentioned in the introduction, our numerical simulations indicate three distinct 
respects in which the evolution of puA differs from the standard qubit evolution given in 



Eq. (61 ). (i) The angular frequency of precession is reduced relative to the standard quantum 
mechanical value Q. (ii) There is a long-term dephasing of the state. For the initial state we 
focus on here, with the spin in the x direction, this dephasing ultimately yields the completely 
mixed state, (iii) The length of the Bloch vector, which would normally maintain a constant 
value of unity, instead oscillates as the spin precesses, achieving its smallest value whenever 
the spin is directed along the y axis. The length of the Bloch vector indicates the purity of 
the state; so we see the state becoming mixed and then regaining its purity every half-cycle 
(long before the purity is reduced permanently because of the dephasing). Evidently what 
is happening is that some of the information in the state is being shared temporarily with 
the environment and then returned to the U A system. When we trace over the environment 
to get puA, this shared component becomes invisible. We call this shared component of the 
state "the ghost part." 

In the following three subsections we present our numerical results for each of these effects. 



B. Reduced precession frequency 

Perhaps the most obvious way in which the precessing spin in our model departs from 
its behavior in standard quantum theory is that the frequency of precession is reduced. 
Numerically we find that the angular frequency depends on the value of our parameter 
A = s/oj, achieving the standard quantum mechanical value Q only as A goes to zero. We 
show an example of the reduced frequency in Fig. [3} which plots the x-component of the 
Bloch vector, b^, as a function of time for a rather large value of A. 
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FIG. 3: The x component of the Bloch vector as a function of time, compared to the standard 
quantum mechanical prediction, for a spin initially in the x direction and precessing around the z 
axis. The dots (red) show the numerical results, and the dashed curve (blue) shows the standard 
quantum mechanical result. In this plot = 200, s = 30 and uj = 100 — we have chosen X = s/co 
to be relatively large to make the effect visible — and we have set = 27r so that time is measured 
in periods of the regular quantum mechanical precession. 

C. Long-term decoherence 

Over a sufficiently long time period, the precessing spin in our simulations decays to a 



stationary mixed state. When the initial state is given by Eq. (62), that is, when the spin 
is initially in the positive x direction, and when the precession is around the z axis, the 
Bloch vector eventually spirals into that axis, so that the final state of the UA system is the 
completely mixed state eg) I a- If instead we take the initial state to have a non-zero 

value of bz, we find that in the evolving state puA{t) the value of bz remains constant, but 
the X and y components of b again spiral into the z axis, so that the UA system finally 
settles into the constant state (l/4)(/c/ ^ Ia + bzlu ® Za)- That is, for any value of bz 
the phase coherence is eventually lost. We find numerically that the decay time depends 
on the environment dimension A^, increasing with increasing before finally approaching 
a constant value when is large. We present an example in Fig. |4| for which the initial 
Bloch vector is 6 = (1/^2,0, 1/^2). We plot there the length b = ^bl + bl + b^ of the 
Bloch vector, which seems to approach the value 1/a/2, consistent with the picture of the 
vector spiraling into the z axis. In making the figure, we chose a value of A^ that shows the 
large-environment limit. 
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FIG. 4: For a spin initially directed at 45° between the x and z axes and precessing around the z 
axis, the length of the Bloch vector decays to 1/^/2 as the vector spirals into the z axis. The dots 
(red) show the numerical results, and the line (blue) is at the level Here N = 1400, s = 30, 

io = 300, and U = 2tt. 

D. The ghost part 

There is also a periodic change in the length of the Bloch vector as a function of time. 
Our numerical results for 6 as a function of time are shown in Fig. |5| in which the spin 
is initially in the positive x direction, and the time axis is in units of the usual quantum 
mechanical precession period 2tx/VL. Notice that the length achieves its minimum value 
twice in each cycle, corresponding to the times when the spin is pointing in the positive or 
negative y direction. Indeed we find that hy never attains the value 1, while remains zero 
as expected. It is as if the Bloch sphere were somewhat flattened along the y axis, so that 
the Bloch vector precesses along the equator of the resulting oblate ellipsoid. 

The shortened Bloch vector associated with the y axis indicates an increased entropy: 
some information has been lost. But it has been lost only temporarily, as it comes back in the 
next quarter-cycle. Again, it appears that some kind of correlation has been temporarily set 
up between the UA system and the environment. While it may be unusual for a correlation 
to automatically reverse itself, one can certainly find instances of such reversal in standard 
physics. While a light pulse is reflecting off a mirror, for example, it is temporarily correlated 
with electrons in the mirror's silver coating, but once the reflection is complete the correlation 
has been undone and the pulse's state is, in the ideal case, as pure as it was before the 
reflection. (The effect we are seeing is similar to the phenomenon of "false decoherence" as 
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FIG. 5: The length of the Bloch vector varies periodically with a frequency equal to twice the 
precession frequency. The dots (red) show the numerical results when the spin is initially in the 
positive X direction. Here = 200, s = 30, uj = 300, and $7 = 2tt. The minima correspond to 
the times when the Bloch vector points in either the positive or negative y direction. Notice that 
with these values of the parameters, the effect is small, with a reduction in length of about half a 
percent. 

described in Refs. [23H26].) 

One might wonder whether under some circumstances the correlation between the UA 
system and the environment might become so thoroughly mixed up within the environment 
that it could never be undone, in which case the entropy of the UA system would have 
increased permanently. This seems to be what happens after many cycles of precession — we 
eventually get the decoherence observed in Subsection IV. C. But can the information get lost 
even when the spin is not precessing? To investigate this question, we have run simulations 
in which, after a quarter-cycle of precession, we turned off the local Stueckelbergian Sua and 
allowed the ubit and environment to continue to evolve according to Seu for thousands of 
precession periods while the UA system remained "frozen" along the y axis. We then turned 
Sua back on and let the system evolve for another quarter-cycle to see whether the Bloch 
vector would regain its full length. Indeed it did — the information had not been permanently 
lost. Similar numerical experiments, in which the precession axis was changed for the final 
part of the evolution, yielded similar results. Now, it is certainly possible to cause the 
correlation to be lost in the environment while the spin is not precessing. It is sufficient 
to change Beu itself in the middle of the numerical run. But normally we can understand 
a time-dependent Hamiltonian as arising from a stationary Hamiltonian acting on a larger 



29 



system. (We include in the system whatever is causing the Hamiltonian to change.) So it 
does not seem unreasonable to assume a time-independent operator Beu ^ we do here, it 
which case it seems that the "ghost part" can be made to return to the U A system, at least 
in the short term. 

E. Anisotropy 

In all of the above simulations, the magnetic field was taken to be along the z axis. We 
associate with this axis the real matrix 1xj®Zai in the sense that this matrix appears in our 
representation of a spin state along the z axis. In the complex representation, we would use 
the usual Pauli matrices X and Z to represent the spin operators along the x and z axes, and 
since these matrices arc real, the xz plane is special in our model in that the associated spin 
matrices in the real representation do not involve Ju. Equivalently, we can say the y axis is 
special, given our conventions, in that it is associated purely with Ju and not with ly. All 
the results we have seen above would be essentially unchanged if we were to have the spin 
precessing around any axis in the xz plane. However, when the precession axis is the y axis, 
all of the above effects disappear. There is neither a periodic nor a long-term change in the 
length of the Bloch vector, and there is no reduction in the frequency of precession relative 
to standard quantum mechanics. For a precession axis intermediate between the xz plane 
and the y axis, the above effects all appear but are not as large as when the precession axis is 
in the xz plane. Consider the decoherence, for example. For an intermediate precession axis, 
we observe that the Bloch vector again spirals into the axis of precession without changing 
its component along that axis. That is, we observe what in the complex theory would be 
called a loss of phase coherence in the energy basis. But the coherence time becomes longer 
as the precession axis becomes more parallel to the y axis. 

There is, however, one effect pertaining to the y axis that has nothing to do with pre- 
cession. Suppose we have no magnetic field — that is, we set Sua to zero — and we choose 
the initial state P(7a(0) to be (l/4)(7[/ ® Ia^ Ju ® Ja)- That is, we try to start the spin 
in the positive y direction. Then one finds that over an extremely short time, the Bloch 
vector shrinks to a shorter length (still in the same direction). This is what one expects 
from Section III: the coefficient of Ju quickly decreases by the factor 1 — A'^. A Bloch vector 
of this length is in fact even shorter than what we would get by starting the spin in the x 
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direction and letting it precess for a quarter cycle. As we will see in the next section, the 
length in the latter case is 1 — In either case, a literal reading of the ubit model would 

seem to say that it is impossible to prepare a pure state of spin in the y direction. Instead, 
we can prepare only a mixed state in that direction. In this sense the Bloch sphere really is 
flattened into an oblate ellipsoid. States that lie beyond the boundary of this ellipsoid are 
simply inaccessible. In Section VI we will introduce an alternative interpretation in which 
a pure state in the y direction is possible, but even in this reinterpretation the predicted 
physics depends on which axis in space we associate with the imaginary Pauli matrix Y . 

Clearly this anisotropy is an effect we need to think about carefully. In the case of 
spin, the special direction — that is, the y direction — is a direction in space. But for other 
physical realizations of a qubit, e.g. a two-level atom, the "direction" one associates with 
the imaginary Pauli matrix Y is not a direction in space; usually it is associated with 
a particular equal-magnitude superposition of the ground and excited states. Moreover, 
we normally have a unitary symmetry that allows us to freely re-express any problem in 
whatever basis we choose — the choice of basis has no physical significance. However, the 
ubit model forces us to treat separately the real and imaginary parts of a Hamiltonian 
or a density matrix — this is why one axis is distinguished as special in the case of the 
precessing spin — and this separation between real and imaginary parts could be changed 
simply by changing the basis. To put it in other words, for any given Hamiltonian there are 
many distinct Stueckelbergians, depending on the basis in which the Hamiltonian is written. 
Section VII explores this question further in the case of higher dimensions. For now, though, 
we try to explain analytically the three effects described in the preceding subsections. 



V. A PRECESSING SPIN ANALYTICAL TREATMENT 



We begin our analysis with Eqs. (47) and (51), in which we have already assumed that 



the state py is continually being projected onto the space of matrices that commute with 
Seu- We write those equations again here: 



dt 



-i[ujKA + iLA,aj]. 



(66) 
(67) 
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For the particular case we are considering now, the local Stueckelbergian is 



Sua = Iu ® La - Ju ® Ka 



-Ju ® Za, 



(68) 



so that Ka = {^/2)Za and La = 0. Inserting this expression into Eq. (67) gives us 



~dt 



■^—{Zaj-ajZ), 



(69) 



where again uj = —{{"^jIIe ® Jul'^j)- We can solve Eq. (69) to get 



j{ujn/2)zt 



(70) 



We are assuming an initial state given by Eq. (62), in which each aj(0) is equal to (1/2)(/a + 
Xa)- In that case we have 

1 



crj{t) = -[Ia + cos{i'jQt)XA + sm{i'jQt)YA] 



(71) 



where Ya = ^ J4 is the usual imaginary Pauli matrix. (Again aj can have a nonzero imaginary 
part. But in p\\ all imaginary contributions will cancel.) The effective density matrix of the 
whole system is 



(72) 



Our strategy will be to try to isolate each of the three effects described above by con- 



sidering different terms of the perturbation expansion of Eq. (72). (i) To see the frequency 



reduction, we expand uj to second order while approximating |\E'j)(\l/j| in Eq. (72) with its 
unperturbed value, (ii) A spread in the values of uj would lead to interference when we do 



the sum in Eq. (72), which would appear as decoherence. But the z/j's begin to diverge from 
each other only in third order. Therefore, to isolate the long-term decoherence, we expand 
Uj to third order while continuing to treat |\I'j)(\l/j| as unperturbed, (iii) The ghost part 
represents a correlation between the ubit and the environment. So to see the ghost part. 



we will expand |\I'j) in Eq. (72) out to lowest nontrivial order (it will be first order) 
while restricting our approximation for uj to second order so as to avoid the complications 
of decoherence. 
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A. Reduced precession frequency 



We begin with Eq. (49) and expand each I'^f) out to second order in A = s/u. Starting 
with the expansion of given in Appendix A, we find that 



-l{^i\lE ® Jf/|vl/±) = ± 1 - - J] |($+ 1^1 >| 



(73) 



where we have replaced the single index j with the pair of indices n and ±. The + and — refer 
to the subspaces in which the operator G takes positive and negative values, respectively. As 
we have done before, for each value of n and k we replace | ($,^ \ V\ ^k)"^ with its ensemble 
average, thereby arriving at 

2 



± 



(74) 



At this order of perturbation theory there is no dependence on n. The factor |\E'j)(\l'j| in 



Eq. (72) we treat as unperturbed, so that ^ can be replaced with ^ |$^) 

P± = {l/2)ilE (S) (/{/ =F Ju)- With these substitutions, one finds that 



1 

AN 



IeIuIa + cos 



1 



A^ 



Qt 



IeIuXa + sin 



1 



A^ 



nt 



IeJuJa 



(75) 



where we have left out the tensor product symbols. Tracing over the environment, we get 
that the density matrix of the U A system is 



PUA{t) 



IuIa + cos 



1 



A^ 



IuXa + sin 



A^ 



JuJa 



(76) 



That is, with these approximations the spin precesses as usual but with its frequency reduced 
by the factor 1 — A^/2. In Fig. [6]we compare this theoretical prediction with the numerically 
observed evolution. 



B. Long-term decoherence 



In Eq. (72), both the EU factor and the A factor depend on j, so that we cannot in 
general separate these two parts of the system when we do the sum. However, to try to 
get an analytic handle on the decoherence, we assume that within each of the two main 
subspaces the EU factor is not significantly correlated with the A factor; so within each 
subspace we can say that the average of the product is the product of the averages. If 

33 



X component 




FIG. 6: Comparison of the numerical data with our analytic result for the x component of the Bloch 
vector, when the spin is initially in the x direction and is precessing around the z axis. The dots 
(red) are the numerical results and the continuous curve (purple) shows our analytic prediction. 
Again the dashed curve (blue) represents the standard quantum mechanical precession. Here 
N = 200, s = 30, w = 100, and Vt = 271. 

we also continue to assume that we can replace each |\E'j)(\l/j| with its unperturbed value 
|<l>j)(%|, we get 



P\\{t) 



4Ar2 



+ P 



^ [I A + cos(z/+^]^)X^ + sin{u+Qt)YA] 

n 

+ cos{iy-nt)XA + sm{u-nt)YA] 

m 

Again using P± = {l/2)ilE ® {lu T Ju), we can rewrite this expression as 

P\\^^^ = ^Y1 [^eIuIa + cos(z/+(]t)/s/c/XA + ?,m{ut^t)lEJuJA] 



(77) 



(78) 



In writing this last equation we have used the fact, mentioned earlier, that the vectors \^ j) 



come in complex-conjugate pairs, and that the real values Vj given by Eq. (49) come in pairs 



with equal magnitudes and opposite signs. This fact is what allows us to combine the two 



sums in Eq. (77) into a single sum. Numerical tests confirm that Eq. (78) yields a very close 



approximation to puA- For example, in Fig. [7]the dashed curve shows the length of the Bloch 



vector as predicted by Eq. (78) — with the values of z/,- determined numerically — while the 



dots represent the numerical results obtained directly from Eq. (64). The good agreement 



provides evidence in support of our assumption of continual projection, made in Subsection 
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III.D, as well as for the specific assumptions leading to Eq. (78) in the present section. Still, 
we would prefer an equation that does not require an exact determination of z/+. So we take 
our approximation a step further. 

In order to evaluate vt to third order, we find it convenient to use the relation 



i{^n\lE ® Ju\^n) =9n- A (^'„ | , 



(79) 



which comes from Eqs. (22) and (23). Thus it is sufficient to expand gn to third order 
and Y^n) to second order. On carrying out this expansion, we find that the third-order 
contribution to — z(\E'^|/s ® Ju\^^) is 

A3 



(•^-^ I y H) vp-v |$+> - ($+1 vp^vp^v \<^i 







\VP-V 






2 


\ n 


v\n)-{^t 


V 





(80) 



We consider the three terms in this expression separately. First, following the reasoning 



in Eq. (29) we approximate ($+| VP-V |$^) as 1, so that the first term square brackets in 



Eq. (80) can be approximated as ($^| V |$^), which we have called f„. We can write the 



second term as 



($+1 VP.VP.V |$+> = 5^ |($+| F 1$, >f {^i\v 1$, 



(81) 



Again we use |($^| |$^)|^ ~ so that we are left with a sum over the eigenvalues 

of the negative-subspace part of V . Those eigenvalues have a typical size that does not 
depend on A^, but their ensemble average is zero, and we expect their sum to be of order 
a/ZV because of random fiuctuations. Thus the whole term diminishes as 1/\/N, and since 
we assume a large environment dimension we take this term to be zero. 



The third term in Eq. (80) is more complicated. We can approximate the numerator 
as The denominator, which we can write as f„ — Vk, can be small, so that the sum 

might depend crucially on the spacing of the values Vk- Those values follow a semicircle 
distribution, but this fact does not tell us how the difference f„ — Vk is distributed. To get 
a somewhat crude approximation, we ignore the issue of the spacing of values and simply 
replace the sum with an integral, assuming a semicircle distribution, and take the Cauchy 
principal value of the integral. This gives us 

l/N _ 1 _1 ... Vr. 



-rj{v)dv 



(82) 
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for the semicircle distribution 77(f) = {N/'K)^Jl — {v/2y. The third term then combines 
quite simply with the first term to give us 

-^{^t\IE^Ju\KY'^-^^■ (83) 



We can now put this result back into Eq. (78) and convert the sum to an integral, again 
assuming a semicircle distribution for f„, obtaining a result for the decay similar to what 
we saw in Section III. 

PUA ~ ^ {IuIa + fit) [cos i^nt) IuXa + sin {^nt) JuJa]} , (84) 



where 

2Ji (A3fit/2) 



fit) 



(85) 



xmt/2 

and ^ = 1 — is the frequency reduction factor we computed in the preceding subsection. 



Fig. 7 shows the length of the Bloch vector as a function of time, as computed from Eq. (84), 



and compares this approximation with the numerical values and with Eq. (78). Clearly our 
approximation of the i/j's is not ideal, but it seems to give us at least a reasonable estimate 
of the coherence time r, which is of order 

^ (86) 



(In fact T has to be over five times this value to make /(r) less than 1/e, though of course 
the curve is not exponential.) Moreover, the detailed shape of the curve traced out by our 
numerical results in Fig. [7] surely depends on the specific model of the environment we have 
chosen, whereas the scaling of the coherence time with A and fl has a better chance of 
carrying over to other models. 



C. The ghost part 



To understand the ghost part, we begin by replacing each uj in Eq. (72) with its second- 
order value, which is ±^ for z/^. (Again ^ = 1 — A^/2.) This approximation allows us to 
write p\\{t) as 



pwit) 



^ ( E K) {K\ ® ^Ht) + E 1^" > i^n I ® ^-(t)] 

\ n n / 



(87) 
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FIG. 7: The spin begins in the x direction and spirals into the z axis as it precesses around 
that axis. Here we plot the length of the Bloch vector as a function of time. The dots (red) 



show the numerical results, and the dashed curve (blue) gives the prediction of Eq. (78) with the 



Uj^s computed numerically. The solid curve (purple) shows our analytic approximation given by 



Eqs. (84) and (85), in which we have made a simple approximation for the values of uj. Here 
N = 1400, s = 30, oj = 300, and = 27r. The "bounce" we see here does not show up in Fig. [4j 
because in that case a small wobble around the z axis does not significantly affect the Bloch vector's 
length. 



where 



= ^ [Ia + cos{Cnt)XA ± sm{^nt)YA] . (88) 



Expanding and to first order in A, we find 



J] |vl>+> (V1/+I = P+ + ^ [P_VP+ + P^VP^] (89) 



and 

A 



J2 \K) {K\=P--^ [P-VP+ + P^VPA . (90) 



Upon inserting these expressions in Eq. (87) we get 



Pii(t) = ^ [IeIuIa + cos{int)lEluXA + sin(efit) {IeJu + A (P-5P+ + P+BP_)} J a] , 

(91) 

where we have again left out the tensor product symbols. 

The term proportional to A is what we are calling the ghost part. Notice that it accom- 
panies what we would normally think of as the y component of the Bloch vector (that is, 
the part proportional to Ju ® J a)- Except for the ghost part and the frequency reduction, 
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the above expression is identical to the standard quantum mechanical evolution given in 



Eq. (61). Of course our expression for the ghost part is valid only to first order in A. If 
we expand each out to second order and trace over the environment, we find that the 
density matrix of the U A system is 



PUA{t) 



IuIa + cos{^nt)IuXA + 1 



(92) 



which shows the shortening of the y component of the Bloch vector. Here we have assumed 
a large environment, so that TteBeu can be taken to be zero. (The ensemble average of 
that partial trace is zero, with fluctuations of order unity. The factor of 1/ (4A^) normalizing 
the density matrix renders such fluctuations negligible.) Under this assumption, the ghost 



part has disappeared. Fig. [8j compares Eq. (92) with our numerical results. 

length of the Bloch vector 
1.000 




time 



FIG. 8: Comparison of our analytic result with the numerical data for the periodic variation of the 
length of the Bloch vector. The spin starts in the positive x direction and precesses around the z 



axis. The dots (red) show the numerical results and the curve (purple) is obtained from Eq. (92). 
Again the minima occur when the Bloch vector is in the positive or negative y direction. Here 
N = 200, s = 30, w = 300, and O = 27r. 

The picture that emerges, then, is that part of the y component of the Bloch vector 
has been lost, but it has been replaced with the ghost part, which represents a correlation 
between the ubit and the environment. When one traces over the environment, what remains 
for the UA system is a Bloch sphere that has been flattened along the y axis by the factor 
1 — A^/2. However, there is an alternative interpretation of the ghost part that seems at 
least as reasonable; this alternative interpretation is the subject of the next section. 
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VI. ALTERNATIVE INTERPRETATION 



Initially we assumed that Alice could perform any measurement on the UA system. One 
such measurement for the case of a spin-1/2 particle is to test whether the spin is in the 
positive y direction. In ordinary quantum mechanics this test would be represented by the 
projection operator (l/2)(/ + Y). The direct translation of this operator into real-vector- 
space terms, according to the prescription of Section II, is the rank-two projection operator 
(l/2)(/[/ ® /yi + J[/ ® Ja)- If this measurement were performed on the completely mixed 
state and the "yes" outcome were obtained, the state of the UA system would be collapsed 
into the state puA = (V4)(-^c/ ^ Ia + Ju ^ Ja), which is the real- vector-space version of spin 
in the positive y direction. But this state cannot persist for any nonzero duration under 
the projection assumption of Subsection III.D. So it would seem that Alice cannot prepare 
a pure state of spin in the y direction, as we noted earlier. 

But there is another way Alice might try preparing a spin state in the y direction. She 
could perform the measurement (1/2) (J^/ ^ Ia + lu ^ ^u) — testing for spin in the positive x 
direction — and upon obtaining the outcome "yes" she could allow the spin to precess around 
the z axis for a quarter-cycle, thereby preparing the effective EU A state 

= ^ [IeIuIa + {IeJu + A {P-BP+ + P+BP^)} Ja] , (93) 



in accordance with Eq. (91 ). Moreover, if at some later time she wanted to test for this state, 
she could do so by first allowing the spin to precess by a quarter-cycle (in the same direction 
as before) and then performing the measurement (1/2)(J{/ ® Ia — lu ® -^a), that is, a test 
corresponding to the negative x direction. (Recall our numerical experiments in which the 
ghost part could be recovered after a long time during which there was no precession.) This 
sequence of operations is perfectly permissible according to our rules; so it should count as 



a valid measurement. In this sense the state py given in Eq. (93) acts like what we would 
normally think of as a pure state of spin in the y direction. One can test for this state and 
get the "yes" outcome with unit probability. So it seems that Alice can prepare a pure spin 
state in the y direction after all. It does not look like a pure state when one traces over the 
environment — it looks like a mixed state — but it acts like a pure state. Again, this effect 
can be seen as an example of "false decoherence" [23H26] , in which part of the environment 
adiabatically follows the evolution of the system of interest. In cases of false decoherence it 
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is misleading simply to trace out the environment, and we seem to have the same kind of 
situation here. 

In our alternative interpretation, then, the Bloch sphere is not flattened. To first order in 
A, a general pure state of spin would be expressed in this interpretation by the EU A density 
matrix 

P^; = ^ [IeIuIa + KIeIuXa + hy {IeJu + A {P-BP+ + P+BP_)} Ja + KIeIuZa] , (94) 



for some unit vector h = {bx,by,bz)- (Eq. (94) is valid only in the special case we have been 



considering in which the environment starts out in the completely mixed state. For a more 
general initial state the form would be different, as we will see in the following paragraphs.) 
Alice can prepare any such state, and if |6| = 1 she can test for the state. It is only that the 
mathematical description of the state is not what we would have expected, since it implicates 
part of the environment. 

We now spell out the alternative interpretation more completely and for a more general 
case. Let A be a d-dimensional system, and let us assume for now that uj can be approx- 



imated by its second-order expansion. (This assumption will be relaxed shortly.) Eq. (47) 
gives the general form of a state consistent with our projection assumption. Our observer 
Alice has no direct control over the environment portion of the state, but according to our 
initial assumptions she can at least prepare a state cr(0) of the system A. If she does so, the 
effective state of the whole system will have the form 

Pll(O) = ^ E^- + ® ^(0), (95) 

n 

where ct(0) is a real, positive semi-definite matrix with unit trace, and the non-negative 
coefficients c„ sum to unity. These coefficients are determined by the initial state of the 
environment and ubit (the initial state of the environment is no longer assumed to be the 
completely mixed state), and we assume that Alice has had no control over their values. 



Alice can manipulate the initial state (95) by choosing a Stueckelbergian Sua of the form 
Sua = Iu ® La~ Ju ® Ka- In general the application of this Stueckelbergian would cause 
the matrix (t(0) to evolve in a different way for each term |\E'^)(\E^^|, but as we have seen, 
to second order in A there are only two distinct value of z/j, namely, = ±C,- So the initial 
state evolves into 

PllW = ^E^- ® a(t) + , (96) 

n 
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where 



(97) 



and the effective Hamiltonian is H' = ^K^ + iL^ in accordance with Eq. (51). 



We can rewrite Eq. (96) as 



a 



a 



) 



(98) 



where 



r^ = E^"l*n)(*;^l (99) 

n 

and we have written a{t) simply as cr. Now, every complex d x d positive semi-definite 



matrix with unit trace can be written in the form given in the right-hand side of Eq. (97) 



for some real a{0) and some Hermitian H'. Because Alice can control a{0) and H' (and t), 
she can determine the matrix a in the state (98). (She cannot control F^.) We can also 



write Eq. ( 98 ) in the following way: 



(X (g) Re (7 + J7 ® Im a) 



where 



X = r+ + r- 



and 



J = ^(Y+- r-) . 



The matrix a evolves according to the equation 

da 

Itt 



[-iH\a]. 



(100) 



(101) 



(102) 



In this respect a behaves like a density matrix. Eq. (100) is reminiscent of Eq. (11) in 



Section II, but the matrices X and J act on the whole EU space rather than just on U . 

Our alternative interpretation is based on an orthogonal transformation O on the EU 
system that turns X and J into tensor products. The intuition is that O identifies a "dressed 
ubit" U' which takes the place of the original ubit in the description of any state of Alice's 
system. Specifically, in Appendix C we construct a transformation O having the following 
effect: 



\n){n\ ® {lu - iJu), 



(103) 



where the vectors \n) constitute an ortho normal (real) basis for the environment. (For our 
purposes it does not matter which basis we choose.) From this equation it follows that 



010 



T 



Pe 



and 



OJO^ = Pe® Ju, 



(104) 
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where pe = ^„c„|n)(n| is a density matrix of the environment. In the alternative inter- 
pretation, our description of Ahce's system is given not by puA but rather by the density 
matrix pu'A) defined as 

Pu>a = T^eOp\\0^. (105) 



When Eq. (100) is vahd, that is, when we can neglect contributions to Vj of higher than 



second order, Eq. ( 104 ) implies that 

1 



Op\\0'^ = ^ pE ® {lu + Ju ® Imcr) 



and therefore 



1 



PU'A 



® Re (T + J(7 ® Im 0") . 



(106) 



(107) 



Note that to this order in perturbation theory, Op\\0^ is equal to the tensor product Pe^Pu'A 
and we can work out the evolution of pu'A without explicitly referring to the environment. 
In fact the theory is almost the same as what we would have gotten simply by applying 
the complex-to-real prescription of Section II. The only difference is the constant factor 
,^ = 1 — in the effective Hamiltonian H' = ^Ka + iLa- 

We now relax the restriction to the second-order expansion of Uj. Starting with the same 



initial state (95), we find that the states Alice can create have the form 



(108) 



where each cr„ is a positive semi-definite complex matrix with unit trace. This equation 



generalizes Eq. (98). Applying our transformation O, we get 



Op\\0^ = ^ '^Cn\n){n\ (/[/ (g) Re (T„ + J[/ ® Ima^) 



(109) 



which is no longer a tensor product between the environment and the UA system. Now 
tracing over the environment gives us 

1 



PuiA = TteOpwO^ = ^{lu 0Rea + Ju Imcr), 



(110) 



where a = Yln'^n^'n- Thus it remains true that Alice can describe her system hj a. d x d 
complex density matrix a. However, in general this matrix will not evolve according to 



Eq. (102). It might decohere, for example. 
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It is interesting to follow the evolution of pu'A numerically for the case of a precessing 
spin. That is, we compute the evolution of p\\{t) as always, but instead of simply tracing 
over the environment to get puA, we apply O before performing the trace, so as to get pu'A- 
As one would expect from the results of this section, the ghost effect is no longer seen: over 
the short term, the Bloch vector associated with pu'A retains its full length over the whole 
cycle. However the other two effects — the reduced precession frequency and the long-term 



decoherence — persist and appear to be the same as before. In Eq. (102) the reduction in 
frequency appears as the constant ^ in H' . We do not see the decoherence in this equation 
because it does not include third-order contributions to Vj. 

The notion of a dressed ubit raises a new question. Throughout our analysis we have 
assumed that Alice is able to build a device that implements an arbitrary Stueckelbergian 
Ie ® Sua (though as we have seen, the ubit-environment interaction can render certain 
components of this operator ineffective). However, if Alice can prepare states of the form 



(100), involving the environment-ubit matrices X and JT, one might wonder whether she also 
has the ability to implement a Stueckelbergian of the form S = X®Sa + J®Ta, where Sa is 
antisymmetric and T4 is symmetric. Evidently she has some control over the environment, 
but does she have the right kind of control to be able to implement SI It seems likely that 



the correct answer is no. Whereas her ability to prepare states of the form (100) follows 
from our initial assumptions, we have not found a way by which she might implement this 
kind of Stueckelbergian. Moreover, the theory is self-consistent as it stands, even though 
Alice's local Stueckelbergians are linear combinations of Ijj and Ju in the original basis, 
while the states she can prepare can be expressed as such linear combinations only after we 
have rotated the EU system by the transformation O. Indeed, Alice prepares these "exotic" 
states by using "ordinary" Stueckelbergians, relying on the interaction matrix Beu to bring 
the environment into the picture. 

Our alternative interpretation also forces us to extend the argument of Subsection III.E 
showing that there can be no signaling through the ubit. If Alice can prepare states involving 
the environment, is it possible that she could leave some mark on the environment that could 
be read by a second observer? The answer is no. If the second observer does not have access 
to the A system itself, then that observer's knowledge is contained in the reduced density 



matrix Tr^py. For the general form given in Eq. (108), this partial trace is always equal to 
(1/2)X, regardless of what Alice has done. So Alice cannot send a signal by this means. 



43 



VII. THE REAL-IMAGINARY SPLIT 



In our analysis of the precessing spin, we began with the Hamiltonian H = h{Q/2)Z and 



converted it into a Stueckelbergian according to the prescription given in Eq. (10), which 
asks us to spht H into its real and imaginary parts. The result was Sua = — (^/2) Jt/ Za- 
But for the very same physical situation, we could just as well have associated the magnetic 
field direction with the Pauli matrix Y instead of Z, in which case the Hamiltonian would 
have been H = h{Q/2)Y and the Stueckelbergian would have been Sjja = (^/2)/(7 Ja- 
This latter Stueckelbergian does not involve the ubit matrix J^, and this fact would have 
made the predictions of our model different from those we obtained in Sections IV and V. In 
particular there would be no decoherence and no retardation of the precession. Moreover we 
would get yet another set of results if we chose to assign to the magnetic field direction some 
other spin matrix such as (X+Y) / a/2. Thus, as we have said before, a given Hamiltonian can 
give rise to many distinct Stueckelbergians, depending on the basis in which the Hamiltonian 
is written. 

Let us see how this ambiguity is expressed for an A system with dimension d. Given a 



certain Hamiltonian, we can imagine changing the basis first and then applying Eq. (10). 
This procedure would give us, for a specified Hamiltonian H, a Stueckelbergian Sua defined 

by 

Sua = /c7 ® Re {-t UHU^) + Ju ® Im {-i UHU^) , (111) 

where ?7 is a unitary transformation representing the change of basis. (After performing this 
operation, we could perform an orthogonal transformation on the A system, but this second 
change of basis does not affect the real-imaginary split.) As the dimension d grows, this 
freedom opens up many new possibilities. For example, even for a set of non-interacting par- 
ticles, the complex-to-real correspondence of Section II could be carried out in an entangled 
basis. And the predicted physics would in general be different for different choices. 

In the specific case of a precessing spin, our model as currently formulated implies that 
in addition to the magnetic field direction, an experimenter should in principle also be able 
to set the orientation of another physically important axis, namely, the axis with which 
we associate the imaginary Pauli matrix Y . The angle between this second axis and the 
magnetic field direction would affect both the rate of precession and the decoherence rate. 

On the other hand, it is conceivable that if something like the ubit model does apply 
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to our world, there could be a law of nature that restricts or determines the way in which 
the real- imaginary split is made. The simplest possibility would be that the conversion 



described by Eq. (10) is required to be carried out in a basis in which the Hamiltonian is 
real. (This law would have meaning only for an isolated system.) Then the Stueckelbergian 
would be uniquely determined by the Hamiltonian — up to an orthogonal transformation on 
the A system — and it would always take the form 

Sua = -Ju ® Ka. (112) 

In this case the ubit would be playing a role very much like that of the phase factor q-^^*/^ 
for a state with definite energy E. In the example of a precessing spin, in our alternative 
interpretation there would be no need to specify a second axis as in the preceding paragraph. 
The reduced precession frequency and the decoherence would always be just as we described 
them in Sections IV and V, because the direction of the magnetic field would always be 
associated with a real spin matrix. Moreover, for an A system of arbitrary dimension, the 
effective Hamiltonian of Section VI would always have the form H' = C,Ka, where Ka is a 
symmetric real matrix — there would be no imaginary part iL^. One consequence is that all 
processes would be slowed by exactly the same factor C,, in which case the slowing would 
be unobservable. Thus if we (i) adopt the alternative interpretation of Section VI and (ii) 



assume that Eq. (10) is always to be applied in a basis in which the Hamiltonian is real, 
then of the deviations from standard quantum theory we have identified, the only one that 
could actually be observed is the decoherence. 



VIII. DISCUSSION 



Over the years a number of researchers have taken an interest in real- vector-space quan- 
tum theory. In a paper published around the same time as Stueckelberg's papers on the 
subject, Dyson argued that when we make use of the time-reversal operator, we are im- 
plicitly basing our theory on the field of real numbers, since that operator is antilinear in 
the complex theory but can be expressed as a linear operator in the real theory [28]. More 
recently Gibbons and others have argued that the complex structure might be an emergent 
property associated with the emergence of a time direction [71 18|. Myrheim has considered 
the possibility of a real-vector-space theory in which Stueckelberg's rule is applied to states 
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but not to observables, and has shown that such a theory allows a canonical commutation 
relation of the form = Jh even for a finite-dimensional state space [20] ■ There has also 
been work on real vector spaces in quantum information theory [TBHT^ , including the 
proof mentioned earlier that n + 1 rebits can simulate n qubits. 

In this paper we have considered a specific model within real-vector-space quantum the- 
ory, characterized by the inclusion of a single binary object, the ubit f/, which is not localized 
and which can interact with anything in the universe. We have focused on characterizing 
the effective theory describing the behavior of a system consisting of the ubit and a local 
real-vector-space object A. From the outset we have assumed that our local observer Alice 
has the ability to prepare any initial state of the UA system and could build an apparatus 
that would implement any Stueckelbergian Sua- However, if both s — the size of a typical 
eigenvalue of ubit-environment interaction — and the ubit's rotation rate u are in effect in- 
finitely large compared to any local frequency, then the local observer's abilities are severely 
constrained. Assuming a large environment, the interaction between the ubit and the envi- 
ronment has the effect of enforcing Stueckelberg's rule on the operators of the UA system. 
If Alice tries to prepare a state puA that does not commute with Jjj^Ia, the noncommuting 
part will instantly disappear. (The noncommuting part has zero trace — so its disappearance 
does not entail a loss of probability.) And if she tries to use a Stueckelbergian Sua that 
does not commute with Ju®Iai its noncommuting part will have no observable effect. Note 
that in our model neither the state p of the whole universe nor the full Stueckelbergian 5* 
commutes with Ie®Ju®Ia- So from a global perspective the theory is quite different from 
standard quantum theory. It is only at the local level that we see all operators commuting 
with Ju- 

Moreover, this automatic enforcement of Stueckelberg's rule at the local level does not 
cause the effective theory to be equivalent to standard quantum mechanics, because the ubit- 
environment interaction has other effects as well. As we have seen, the degree of divergence 
between the ubit model and standard quantum theory depends on how one interprets the 
ubit model. We have distinguished two interpretations, which we have called the literal 
interpretation and the alternative interpretation. 

The literal interpretation assumes that puA is the full and correct description of the state 
observed by Alice. If we adopt this interpretation, then in addition to losing the states and 
transformations that would have violated Stueckelberg's rule, we also lose certain states and 
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transformations that quantum mechanics normally allows. For example, for the spin of a 
spin- 1/2 particle, the Bloch sphere is in effect flattened into an oblate eUipsoid. Except in 
a certain preferred plane, it is impossible to prepare a pure state of spin in any direction. 
In particular, in standard quantum theory there will always be one direction with which 
we associate a purely imaginary spin matrix (usually called the y direction). This is the 
direction in which, according to the literal interpretation, the purity of a spin state is the 
most limited. Moreover, if one prepares a pure state in the favored plane and allows the 
spin to precess out of that plane, then the purity of the state will temporarily be diminished 
in order for the Bloch vector to be able to fit inside the flattened sphere. 

One can sec how the ubit model thus provides a kind of interpolation between standard 
complex-vector-space quantum theory and standard real- vector-space quantum theory with 
no ubit. In the former case, we have the full Bloch sphere. In the latter case, we have 
a "Bloch disk," as all rebit states are confined to a plane. In the ubit model, in the ht- 
eral interpretation, the width of the narrow axis of the flattened Bloch sphere depends on 
the parameter A = s/ui. As this parameter gets larger — e.g., as the ubit's rotation rate 
diminishes — the Bloch ellipsoid becomes more flattened. Indeed, when the ubit is not ro- 
tating at all, there is nothing preventing the ubit's interaction with the environment from 
always bringing puA to the form {l/2)Iu®pA-i in which case the ubit becomes inconsequential 
and, since pa is real, the spin is confined to the Bloch disk. 

However, there are good reasons to prefer the alternative interpretation. In this inter- 
pretation the matrix puA does not describe what Ahce sees. Instead, the state of Alice's 
system is described by pwA = Tr OpyO^, where O is an orthogonal transformation on the 
EU system. In the case of a spin-1/2 particle, in the alternative interpretation the Bloch 
sphere is not flattened: the set of possible density matrices pwA includes all the pure states 
on the sphere. In effect the density matrix pijiA includes the ghost part that would have 
disappeared upon simply tracing over the environment to get puA- It makes sense to de- 
scribe Ahce's states in terms of pu'A, because this description reflects what Ahce is actually 
able to do — what states she can prepare and test for — according to our original assumptions. 
We have noted that though Alice's actions involve the environment, they do so in a way 
that does not convey information. 

The alternative interpretation yields a greater degree of agreement between the ubit 
model and standard quantum theory. To second order in our parameter A = s/ou, the 
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only difference between the two theories, with regard to the behavior of the local system, 
is a factor ^ < 1 multiplying the real part of the Hamiltonian. Moreover, we have also 
contemplated the possibility of a law of nature that would force the local Stueckelbergian to 
be of the form —Ju® Ka, in which case the factor ^ becomes a universal retardation factor 
and therefore unobservable. 

As we go to higher order in A the correspondence breaks down further, and the ubit 
model will surely be very different from standard quantum theory once A is of order unity. 
One difference we have been able to identify is the spontaneous decoherence, which appears 
as a third-order effect. For the case of a precessing qubit, we have an analytic estimate 
of the coherence time in our model: it is r ^ l/(A^i7), where Q is the angular frequency 
of precession. Thus experimental observations of the product tQ can be used to place an 
upper bound on A. A particularly high value of tQ was obtained recently by Chou et 
al. for an optical atomic transition [30]. These authors observed a coherence time around 
10 seconds for a transition with angular frequency = 7 x 10^^ s^^. This result implies that 
our parameter A cannot be much larger than 10~^. Note that the slowing by the factor ^, 
being a second-order effect, could conceivably place a tighter bound on A, but because of our 
uncertainty about the real- imaginary split, it is hard to know how to look for this slowing. 

Again, the model we have described here is highly nonlocal, and it is interesting to ask 
whether we would get similar results from a local theory in which the ubit is replaced with 
a ubit field. (One is reminded of Refs. [IHl HH], in which the authors show that a complex- 
vector-space quantum computation can be simulated locally by a real- vector-space quantum 
computation.) We have relied on the environment to effectively cut out the parts of the 
local Stueckelbergian and local density matrix that do not satisfy Stueckelberg's rule. In 
order to get this kind of environmental effect with a ubit field, we would probably need the 
field to be such that a change at one location would quickly cause similarly large changes in 
the surrounding region, so as to bring a large part of the environment into the dynamics to 
create effects similar to what we have seen here. 

Clearly the work presented here leaves many questions unanswered. It is only an initial 
exploration of the model. But we can at least see that the model is not in obvious conflict 
with observation as long as A is sufficiently small. Note that in a complete theory it should 
not be necessary that this parameter remain constant for all time. So the degree of deviation 
from standard quantum theory could conceivably be different at different stages of the 
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universe's evolution. The effective theory might be experimentally indistinguishable from 
standard quantum theory at the present time but could look much more like ordinary real- 
vector-space quantum theory (with no ubit) at some very early time. In this sense the 
familiar complex structure of quantum theory could indeed be an emergent feature in the 
ubit model. 
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APPENDIX A: BASIC PERTURBATION EXPANSION 

To summarize our notation: the Hermitian operator on which we base our perturbation 
expansions is 

G = Go + AV, (113) 

where X = s/u and the unperturbed matrix Go has only two distinct eigenvalues, ±1, each 
corresponding to an A^-dimensional subspace. The states |$^) are the eigenstates of Go 
that diagonalize V within each of these two subspaces, and are the corresponding 

eigenstates of G. The exact eigenvalues of G are labeled ±gn- 

To do our perturbation calculations, we need to expand gn out to third order and 
out to second order. We write here the expansions for Qn and 1^1/^). As we have mentioned 
earlier, the eigenvector I^P^) is the complex conjugate of I^E'^). Our expansion for the 
eigenvectors is derived from Eqs. (141-143) of Ref. [31]. Our expansion for the eigenvalues 
is obtained from the expanded eigenvectors together with Eq. (7) of the same paper. The 
normalization convention in that paper is such that the expanded vectors are not of unit 
length, but in the expansion of I^E'^) below we have renormalized the vector to unit length 
(to second order in A), which introduces the two second-order terms proportional to |$^). 
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Here we use the notation Mnk as an abbreviation for ($+|M|$^), for any matrix M. Again, 
P_ is the projection onto the space spanned by the states |$~). 
The expansion for gn is 



gn^l + XVnn + Y(VP-V)r 



+ 



A3 



(114) 



Vnn — Vkk 

where the prime on the summation sign indicates that k runs over all values other than n. 
The expansion for |^+) is 

' 2 [ ^'^^ Vnn-V,, 



I*" 



+ 



p_yp_y|$+) - VnnP-vm 



^iVP^V)nn\^:) 
{VP_V)UVP-V) 



+y;(^^^^^p-v 

, V Vnn Vkk 
k ^ 

l \{VP.V)knf 
2 {Vnn — VkkY 

{VP-V)knVnn 
Vnn — Vkk 



1$:) + 



(Vnn — VkkY 
{VP.VP.V)kn 



(115) 
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' {VP.V)kl{VP.V)in 
^ {Vnn-Vkk){Vnn-Vu) 



APPENDIX B: LETTING u APPROACH INFINITY 

For most of this paper we have assumed that both s and io are very large compared to 
Sua- We have also assumed that the environment dimension N becomes arbitrarily large. 
In this Appendix we consider a different limit. Here s and N both remain finite, and we let 
lo go to infinity. Our aim is to show that in this limit, the evolution operator in the ubit 
model, c^*, becomes equivalent to a corresponding evolution operator of standard quantum 
theory. 

Again, our Stueckelbergian for the whole EJJ A system is 

S = -ujIe ®Ju®Ia + sBeu ®Ia + Ie® Sua- (116) 
Let us now write this operator as 

S^-ojJ + b, (117) 
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where J = Ie ® Ju ® Ia and D = sBeu ® Ia + Ie® Sua- 

We begin by rewriting e'^* as follows, using the fact that S is antisymmetric so that 



sm 



= S ^ ^ +cos(Vg^gt). 

VS^S 

Note that ^/ S'^S is well defined since S'^S is a positive semi-definite matrix. Now we make 
the substitution S = —uJ + D: 

s^s = ( - ioj^ + b^){-ujj + b) = io^i + uj{jb + bj) - b^. (119) 

Here we have used the fact that both J and D are antisymmetric. Recall that the part of 
b that commutes with J can be written as be = (1/2) (£) — jbj^ . In terms of be, we have 

S^S = uj^I + 2ujjb, - (120) 

so that 



'121] 



We can therefore write the evolution operator (Eq. (118)) as 



sin(a;ty^/+^-g. / , ^JD, 



e^* = - J + - \ — ' + cos ujt\ I + - ^ . (122) 



When u is very large, we can expand the square roots in powers of 1/u. Ignoring terms 
of order l/u multiplying the sine, and also ignoring terms of order 1/u inside the sine and 
cosine, we get 

e^* ^ - J sin [{ul + Jb,)t] + cos [{ool + jb,)t] = e^-J'-+^^^' = e--^"*e^=*. (123) 

Thus, when u is very large, the dynamics effectively separates into two parts: (i) the ubit 
rotates very rapidly, and (ii) the whole system evolves according to the Stueckelbergian Dc- 
(These two operations commute with each other.) 
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We now want to show that this evolution amounts to an ordinary quantum mechanical 
evolution. The Stueckelbergian Dc is exactly what one would get by starting with the 
Hamiltonian 

^=(^c)oo + ^(^c),o (124) 

and simply rewriting the same physics in real-vector-space terms as in Section II. Here the 
indices and 1 are ubit indices as in that section. When we write out the definition of 
this Hamiltonian becomes 

H ^Hs®Ij^ + Is® H_A, (125) 

where Hs = s{ [{Beu)c\qo + ^ [(^eu)c\ iq} ^nd Hj^ = [('^c/a)c]oo + ^ [('^t/A)c] we are 

using script letters to refer to systems described in terms of a complex vector space. (Again 
the subscript c means that we are taking only the part of the operator that commutes with 
J.) Thus the environment and the local system evolve independently, each according to its 
own Hamiltonian. 

We conclude that, even without letting s or N go to infinity, dynamics in the ubit model 
reduces to the dynamics of standard quantum mechanics as ou approaches infinity. 

APPENDIX C: CONSTRUCTING THE ORTHOGONAL TRANSFORMA- 
TION O 

In this Appendix our aim is to find an orthogonal transformation O on the EU system 
such that 0|\E'^) (\E'^|0^ = (l/2)|n)(n| {lu ~ 'iJu)^ where the real vectors |n) constitute 
an orthonormal basis for the environment. We construct O in two steps. First, let Oi be 
the matrix 

Oi = E(l^n)(*:i + l^;)(*;l)- (126) 

n 

Oi transforms between two orthonormal bases, so it is unitary. It is also real and therefore 
orthogonal. Upon applying Oi to |^'+)(^+|, we get 

= |$+)($+| = |+)(+| = 1 ® {lu - iJul (127) 

where we are using the factorization |$+) = ® \+) introduced in Subsection III.C. 
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Next, let U be the unitary matrix that takes to \n). That is, U = J2n \^){^n\- -^^^ 
let O2 be the real-vector-space version of U according to the transcription rules of Section 
II. That is, 

02 = ReU ^ lu + lmU ^ Ju- (128) 
By rewriting f^l*/*^) (</'^|f^^ = |'^)('^| in real-vector-space terms, we get 

O2 (Re 10+) (0+1 ®Iu + lm 10+) (0+1 ® Ju) Ol = \n){n\ ® lu- (129) 

Multiplying both sides of this equation by —He ® Ju (which commutes with O2) gives us 



O2 (ilm| 



2 Re 



Ju) O2 = —i \'n){n\ (g) Ju. 



Now we add Eqs. (129) and (130) to get 



O2 [\<Pt){<Pt\ ® (lu - ^Ju)] O2 = \n){n\ ® - iJu) 



(130) 



(131) 



in which the left-hand side mirrors the right-hand side of Eq. (127). Finally we define O to 



be O = O2O1. Then Eqs. (127) and (131) imply that 



0\^l){n\0^ = i^\r>){n\®{Iu -iJu). 

which is what we wanted to show. 

It is interesting to work out the effect of O on Seu- We can write Seu as 



(132) 



Seu = -i^Y.3n (l^^)«l " l*;)(^;l) , (133) 

n 

where the QnS are eigenvalues of G as in Subsection III.B. Upon applying O, we get 

OSeuO'^ = -uj(^Y^gn\n){n\^ Ju. (134) 



Thus the transformation O brings Seu to a tensor-product form. Any antisymmetric real 
matrix such as Seu can be brought to block-diagonal form by an orthogonal transformation, 
with 2x2 blocks proportional to J [27]. If we take the basis defined by the |n)'s to be the 
standard basis for the environment, then O is an orthogonal transformation that takes Seu 
to this form. 



[1] G. Birkhoff and J. von Neumann, Ann. Math. 37, 823 (1936). 



53 



[2] D. Bohm, Quantum Theory (Prentice-Hall, New York, 1951). 

[3] E. C. G. Stueckelberg, Helv. Phys. Acta 32, 254 (1959). 

[4] E. C. G. Stueckelberg, Helv. Phys. Acta 33, 727 (1960). 

[5] A. Trautman, in On Einstein's Path (Springer, New York, 1966), p. 487. 

[6] P. J. Lahti and M. J. Maczynski, J. Math. Phys. 28, 1764 (1987). 

[7] G. W. Gibbons and H. J. Pohle, Nucl. Phys. B 410, 117 (1993). 



[8] H. A. Chamblin and G. W. Gibbons, |ar xiv:gr-qc/9510006 (1995). 
[9] J. B. Barbour, Phys. Rev. D 47, 5422 (1993). 



[10] L. Hardy, |a^xiv:quant-ph/0101012| (2001). 

[11] C. M. Caves, C. A. Fuchs, and R. Schack, J. Math. Phys. 43, 4537 (2002). 
[12] S. Aaronson, |arxiv:quant-ph/0401062| (2004). 
[13] P. Goyal, New J. Phys. 12, 023012 (2010). 

[14] G. Chiribella, G. M. D'Ariano, and P. Perinotti, Phys. Rev. A 84, 012311 (2011). 
[15] S. Kobayashi and K. Nomizu, Foundations of Differential Geometry (Interscience, New York, 
1963-1969), V. 2, Chap. IX, Sec. 1. 



[16] T. Rudolph and L. Grover, arxiv:quant-ph/0210187 (2002). 



[17] J. M. Fernandez and W. A. Schneeberger, |arxiv:quant-ph /0307017 (2003). 
[18] M. McKague, M. Mosca, and N. Gisin, Phys. Rev. Lett. 102, 020505 (2009). 



[19] M. McKague, arXiv:1109.0795| (2011). 



[20] J. Myrheim, [arxiv:quant- pti79905037 (1999). 

[21] S. D. Bartlett, T. Rodolph, and R. W. Spekkens, Rev. Mod. Phys. 79, 555 (2007), Eqs. (4.25, 
4.26). 

[22] E. Lubkin, J. Math. Phys. 19, 1028 (1978). 

[23] A. J. Leggett, in Applications of Statistical and Field Theory Methods to Condensed Matter 
(Proceedings of the 1989 NATO Summer School, Evora, Portugal), edited by D. Baeriswyl, 
A. R. Bishop, and J. Carmelo (Plenum Press, 1990). 

[24] A. J. Leggett, in Quantum Measurement: Beyond Paradox, edited by R. A. Healey and G. Hell- 
man (U. Minnesota Press, Minneapolis, 1998). 

[25] J. R. Anglin, J. P. Paz, and W. H. Zurek, Phys. Rev. A 55, 4041 (1997). 

[26] W. G. Unruh, in Relativistic Quantum Mechanics and Decoherence, edited by H. P. Breuer 
and F. Petruccione (Springer, New York, 2000), p. 125. 



54 



[27] R. A. Horn and C. R. Johnson, Matrix Analysis (Cambridge Univ. Press, 1985), p. 107. 

[28] F. J. Dyson, J. Math. Phys. 3, 1199 (1962). 

[29] L. Hardy and W. K. Wootters, Found. Phys. 42, 454 (2012). 

[30] C. W. Chou, D. B. Hume, M. J. Thorpe, D. J. Wineland, and T. Rosenband, 

Phys. Rev. Lett. 106, 160801 (2011). 

[31] J. O. Hirschfelder and P. R. Certain, J. Chem. Phys. 60, 1118 (1974). 

[32] A linear operator with this property is said to define a complex structure on a real vector space 

mm- 

[33] Alternatively, we could write p = (1/2) [(1/2) (/t/ - Yu) a + (1/2) (It/ + Yu) ® a], where Y 
is the imaginary Pauli matrix and the bar indicates complex conjugation. This form is an 
extension of the general prescription given in Ref. |21) for bringing a reference frame into the 
quantum description. In this case the ubit would be the reference frame. We do not use this 
formulation here because we want to emphasize the possibility of treating the ubit as a rebit 
rather than as a qubit. For a rebit, we cannot have (/ it Y)/2 as possible states along with 
the states (/ it X)/2 and (/ it Z)/2, which we will need. 

[34] We say two matrices Mi and M2 are perpendicular when Tr Af|M2 = 0. 



55 



